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Context. Gamma-ray bursts (GRB) are powerful, short duration events with a spectral luminosity peaking in the keV - MeV (BATSE) 
range . The prompt emission is thought to arise from electrons accelerated in internal shocks propagating within a highly relativistic 
outflow. 

Aims. The launch of Fermi offers the prospect of observations with unprecedented sensitivity in high-energy (HE, > 100 MeV) gamma- 
rays. The aim is to explore the predictions for HE emission from internal shocks, taking into account both dynamical and radiative 
(~| \ aspects, and to deduce how HE observations constrain the properties of the relativistic outflow. 

^ Methods. The prompt GRB emission is modelled by combining a time-dependent radiative code, solving for the electron and photon 

I ' distributions, with a dynamical code giving the evolution of the physical conditions in the shocked regions of the outflow. Synthetic 

O ' lightcurves and spectra are generated and compared to observations. 

Results. The HE emission deviates significantly from analytical estimates, which tend to overpredict the IC component, when the time 
C/3 , dependence and full cross-sections are included. The exploration of the parameter space favors the case where the dominant process 
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in the BATSE range is synchrotron emission. The HE component becomes stronger for weaker magnetic fields. The HE lightcurve can 
display a prolonged pulse duration due to IC emission, or even a delayed peak compared to the BATSE range. Alternatively, having 
dominant IC emission in the BATSE range requires most electrons to be accelerated into a steep power-law distribution and implies 
strong second order IC scattering. In this case, the BATSE and HE lightcurves are very similar. 

Conclusions. The combined dynamical and radiative approach allows a firm appraisal of GRB HE prompt emission. A diagnostic 
procedure is presented to identify from observations the dominant emission process and derive constrains on the bulk Lorentz factor, 
particle density and magnetic field of the outflow. 
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^ ■ 

I 1. Introduction high-energy spectral indices and break energies of GRBs that 

^— ^ ' have significant MeV emission. More recently iGiuliani et al.l 

• • ■ The forthcoming first results of the Fermi gamma-ray space ,2008) reported observations of GRB 0805 14B by AGILE 

, fe/e^cope call for a detailed study of the high energy (above 100 showing some evidence that the emission above 30 MeV 

MeV) gamma-ray burst (GRB) emission. Current observational ^^^^^^^ ^ j^^g^^ j^^^^i^^ ^^an the emission observed at 

H , information on very high-energy gamma-rays emitted in a GRB j^^g^ energies. Evidence of even higher (TeV) energy emission 

,date from the EGRET (Energetic Gamma-Ray Experiment from GRBs was reported from ground-based experiments, based 

Telescope) mission on board the CGRO (Compton Gamma Ray detection of extensive air showers pr oduced by high en - 

Observatory). It detected high energy photons from a handful of g^gy photons propagating in the atmosphere (lAtkins et al.l2000l) . 
GRBs; the most energetic (18 GeV) phot on was detected in the 
case of GRB 940217 (Hurley et al, 1994). , Gonzalez et al.l(l2003h 

reported the observation of a bright high-energy component The observation of high energy spectral components in 

in GRB 941017, showing a strong temporal evolution, distinct GRBs can provide strong constraints on present models for 

from the low-energy (<2 MeV) component. The inspection of the GRB prompt phase. GRBs are believed to be produced 

the sample of gamma-ray bursts that were observed both by by ultra-relativistic (F ^ 100) outflows ejected from a newly 

EGRET and BATSE (Burst and Transient Source Experiment) formed compact stellar mass source. The prompt gamma-ray 

indicates that these burst s were am ong the brightest ones emission is usually interpreted as radiation from accelerated 

detected by BATSE (e.g. Baring"2006); as BATSE trigger was electrons in shock wave s that propagate within the outflow 

sensitive in the lower energy range, there could be a population (iRees & Meszaros I |1994|) . Such internal shocks can form if 

of bursts with high energy photons that did not trigger BATSE the ejection process by the central source is highly vari- 



dJones et al.l Il996l) . iKaneko et al.l (|2008|) reported the spectral able. A high energy spectral component is expected within 
analysis of combined BATSE and EGRET data for 15 bright this framew ork ( s ee e.g. Papathanassio u & MeszarosI 19961 ; 
GRBs in energy range -30 keV-200 MeV, emphasizing the im- Sari & Piran 1997; Pifla & Loeb 1998; Guetta & Granot 2003 



portance of such broadband spectral analysis in constraining the Pe'er & Waxman 2004; Razzaque et al. 2004; Asano & Inqu 

2007; Gupta & Zhang 2007; iGalh & G uetta 2008; Fan & Pira 

Send offprint requests to: F. Daigne (daigne@iap . fr) 120081: lAndo et al.»2008) . The typical GRB spectrum in the low 
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gamma-ray range, as observed for instance by BATSE, is a 
smoothly connected broken power law with a break energy in 
the range 0.1-1 Me V. This component can be directly produced 
by synchrotron radiation from the shock accelerated electrons, 
or by inverse Compton scatterings of low-energy synchrotron 
photons by the relativistic electrons. Thus observations of 
the GRB spectrum extending to very high energy emission 
(GeV ranges) can be expected when the keV-MeV photons 
are inverse Compton scattered (provided that the jj opacity 
in the source is low). Depending on the relevant parameters, 
the flux of the high energy component can be even compara- 
ble to the prompt GRB gamma-ray flux in BATSE energy range. 



Significant observational progress is expected with the 
launch of Fermi, whose two instruments, GBM {GLAST burst 
monitor) and LAT (Large Array Telescope) will allow the 
observation of GRBs in an unprecedente d spectral range from 
8 keV to 10 GeV or above (.Gehrels & Michelsonlil99 9). The 
LAT has a large field of view (~ 2 sr), is about 10 times 
more sensitive than EGRET and has a very short dead time 
of ~ 100/is (compared to 100 ms for EGRET). Fermi should 
therefore detect 100 to 200 GRBs per year (GBM+LAT), with 
an appreciable number of them being bright enough above 100 
MeV to allow a good characterization of their temporal and 
spectral properties in the high-energy gamma-ray range. 



2. Internal shocks: dynamics and radiative 
processes 

We assume that a relativistic outflow is ejected by the cen- 
tral source of the gamma-ray burst, and that, due to initial 
variability in the distribution of the Lorentz factor, shock 
waves form and p ropagate within this outflow (internal shocks, 
iRees & Meszaros|[T994,) . A fraction of the kinetic energy which 
is dissipated in the shock waves is radiated and produces the 
observed prompt GRB. Here, we focus on the most discussed 
version of the internal shock model, where the radiation is due 
to shock-accelerated electrons in optically thin conditions. It 
has been suggested that shock accelerat ed protons could also 
contribute to the high-energy emission dRazzaque et al.l 120041: 
lAsano & Inoud 120071: lAsano et aLlbOOSl) . or that the emission 
could occur in o ptically thick regions leading to quasi-th ermal 
comp t onization ( Ghisellini & Celotti 1999; Meszaros & ReesI 
I2OOOI: iPe'er & WaxmanI 12004: .Giannios & Spruili ,20071) . or 
that the dominant process is not relat ed to shock-accelerate d 
electrons but rather to decaying pions (Pa czvnski & Xul [19941) . 
These alternative possibilities are not considered in this paper. 

In order to follow the time evolution of the photon spectrum 
emerging from the relativistic outflow during the internal shock 
phase, several steps are needed: 



This paper presents a detailed calculation of the GRB prompt 
emission in the context of the internal shock model, focussing 
on the high energy (above 100 MeV) range. The emission in 
the shocked region is computed using a radiative code that was 
developed to solve simultaneously the time evolution of the elec- 
tron and photon distribution, which is a significant improve- 
ment compared to studies based on an analytical estimate of 
the spectrum. Th is radiative calculation i s less detailed than in 
previous studies ( Pe'er & Waxmanll2005l: lAsano & Inouell2007h 
as it does not include components emitted by shock-accelerated 
protons or by electron-positron pairs created from yy annihila- 
tion. However it includes all the relevant processes for the emis- 
sion from shock-accelerated electrons, whose contribution is ex- 
pected to be dominant. In addition, this radiative calculation is 
combined for the first time with a detailed dynamical simula- 
tion of the internal shock phase, which allows us not only to 
estimate the spectrum of the prompt GRB emission, but to gen- 
erate full synthetic GRBs with lightcurves and spectra. This ap- 
proach is described in Sect.|2] The effect of the parameters de- 
scribing the physical conditions in the shocked medium on the 
shape of the emitted spectrum in the comoving frame is shown 
in Sect.|3] The parameter space of the internal shock model is 
explored in a systematic way in Sect.lH which allows the iden- 
tification of the different classes of high-energy spectra that can 
be expected. We show how Fermi data will allow us to diagnose 
the dominant radiative process (synchrotron radiation vs inverse 
Compton scatterings), the physical conditions in the shocked 
medium (electron distribution, magnetic field) and the proper- 
ties of the relativistic outflow (Lorentz factor and injected kinetic 
power). Finally, Sect. |5] describes examples of synthetic bursts 
(lightcurves and spectra) and discuss how the comparison be- 
tween the LAT and the GBM lightcurves and the observed spec- 
tral evolution in Fermi bursts are also powerful tools to better 
constrain the physical processes at work in GRBs. Sect.|6] sum- 
marizes the results of this study. 



1 . The dynamics of the internal shock phase must be followed 
to determine the physical conditions behind each shock 
wave; 

2. In the shocked medium, electrons are accelerated and the 
magnetic field is amplified. The emitted photon spectrum has 
to be computed from the time-dependent evolution of the rel- 
ativistic electrons. This evolution is governed by several ra- 
diative processes that are in competition with the adiabatic 
cooling due to the spherical expansion. 

3. From the evolution of the emission in the comoving frame of 
the shocked material, one can deduce the observed prompt 
GRB lightcurve and spectrum. 



Some aspects of this project have already been studied by sev- 
eral authors, who focus on the second step (radiative processes 
in the comoving frame) after assuming a typical collision be- 
tween two relativistic shells for the first step. This has been done 
either using an approximate analytical or semi-analytical esti 
mate of the spectrum ( e.g. Papathanassiou & Meszaros 1996 
Guetta & Gra"no3 l2003t l&ipta & Zhang 2007; Galli & GuetL 
20081; iFan & Piranll2008l; lAndo et al.ll2008h or a detailed radia- 



tive code dPe'er & Waxman I2004t lAsano & Inouell2007l) . Such 
studies allow to discuss dififerent emission mechanisms of high 
energy photon production during internal shocks and to de- 
rive the expected high energy photon spectrum from one sin- 
gle shocked relativistic shell. However, they cannot produce full 
lightcurves and time-evolving spectra, and evaluate the role of 
the dynamics of internal shocks in the observed spectral evolu- 
tion. In this work we attempt to improve this approach by com- 
bining a complete model for the dynamics of the internal shocks 
with a detailed calculation of the relevant radiative processes oc- 
curring in the shocked medium. This allows us for the first time 
to obtain the time evolution of the high-energy gamma-ray emis- 
sion in a GRB. The procedure we have adopted in described in 
the present section. 



Bosnjak, Daigne & Dubus: Prompt HE emission from GRBs in the internal shock model 



3 




0.2 0.4 0.6 0.8 1 1 2 3 

M / M,„,,, t„,, / (1+z) [s] 

Fig. 1. Dynamics of internal shocks: an example. Left. Evolution of the distribution of the Lorentz factor in the relativistic outflow. 
In this example, the initial distribution (thick solid line) corresponds to a case where material has been ejected for = 2 s by the 
central source, with a Lorentz factor r(f) increasing from 100 to 400. Two internal shock waves form when the faster part catches up 
with the slower one, as it is shown by the evolution of the distribution of the Lorentz factor with time (thin solid lines). The dynamics 
is computed assuming an injected kinetic power E - 5 y. 10^^ erg s"'. Right. Corresponding physical conditions (T,, p, and e,) in 
the shocked material (see text), plotted as a function of the photon arrival time. For each curve, the two branches correspond to 
the two shock waves that form in the outflow. Two additional quantities are also plotted, Fm and B' , assuming - 1/3, f = 10"^, 
p - 2.5 and - 1/3. Circles indicate the physical conditions at times f = 6.9 x 10^, 1.9 x 10"*, 5.2 x 10"* and 1.8 x 10^ s, for which 
the Lorentz factor distribution in the outflow is plotted in the left panel. Notice that the shock propagating forward reaches rapidly 
the front edge so that, in this example, the emission will be dominated by one shock wave only. This example will lead to a single 
pulse burst. In practice, the initial Lorentz factor considered here has to be seen as a building block for more realistic distributions 
leading to multi-pulse Ughtcurves (see Sect. |5.2t . 



2. 1 . Dynamical evolution during the internal shock phase 



The dynamics of i nternal shocks within a relativistic outflow has 
been described in'Kobavas hi et al.1 (1 1 997') in the case where the 
central engine is emitting a discrete number of shells, separated 
by short periods without any ejection. In this scenario, each 
pulse observed in the GRB lightcurve is due to a collision 
between two shells. One potential problem with this approach 
is that the pulse shape in the decay phase is dominated by 
the so-called curvature effect, i.e. the spreading of the arrival 
time of photons emitted simultaneously on a curved surface. 
Such a decay is too fast co mpared to observations (see e.g. 
ISoderberg & Fenimord l200l]) . In this paper, the dynamics of 
internal shocks is rather computed using the model developed by 
iDaigne & MochkovitchI ([1998.) . where the relativistic ejection 
is now considered as a continuous process. Instead of colli- 
sions between discrete shells, internal shocks are in this case 
shock waves propagating within the outflow. In the observed 
lightcurve, the shape of pulses in their decay phase is then 
determined by the hydrodynamical timescale associated with 
the propagation of the shock waves, rather than the curvature 
effect (except at the very end of this dynamical phase). Slow 
pulse decays can easily be obta ined, which greatly improve s 
the agreement with observations dPaigne & Mochkovitchl2003l) . 



The dynamics during the internal shock phase is entirely 
determined from the following parameters : the total duration 
t^f, of the relativistic ejection and the history of the Lorentz 
factor r(f) and of the injected kinetic power E(t) during this 
ejection. In practice, the outflow is described as a series of 
shells emitted regularly over a timescale At <K so that the 
number of shells is much larger that the number of pulses in the 
lightcurve. These shells interact only by direct collisions, so that 
the propagation of a shock wave is discretized by a succession 
of shocks between shells. The details of the implementat ion of 
this model are described in lDaigne & MochkovitchI (Il998h . This 
method has been validated by a comparison with the results of a 
ID Lagrangian relativistic hydrocode (Daigne & Mochkovitc^ 
2000). Relativistic hydrodynamical sim ulations of i nternal 
shocks have also been pe rformed bv,Mimica et alj (|2004|) in the 
context of blazars and bv lMimica et al. ( 20071) in the context of 
GRBs. The authors discuss the efficiency of the conversion of 
kinetic energy into radiation, and especially the impact of the 
possible magnetization of the outflow, which is not considered 
in the present paper The output of a simulation of the internal 
shock dynamics is the time evolution of the physical conditions 
in the shocked medium behind each shock wave (comoving 
mass density p», comoving specific energy density e», and 
Lorentz factor F*). This is illustrated in a simple example shown 
in Fig.[Tl where the Lorentz factor distribution in the outflow is 
plotted at different times f, and the physical conditions in the 
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shocked medium are plotted as a function of fobs /(I +z) - t-R/c 
(arrival time in the observer frame of photons emitted along the 
line of sight at radius R and time t). 

To estimate the typical radius and shock conditions 
in internal sho cks, a simple "two sh e lls" model is often 
used (see e. g. Rees & Meszaros 1994; Barraud et al. 
iDaigne & Moc hkovitch 2007; Kumar & McMahon 2008). We 
consider the ejection of two equal mass relativistic shells with 
Lorentz factor Ti and r2 from the central source. Shell 1 is 
ejected first and shell 2 after shell 1, with a delay r. If the con- 
trast K - Fa/Fi is greater than unity, an internal shock will occur 
at a radius 



8k^ 



PCT. 



(1) 



where the average Lorentz factor is f = (Fi + F2)/2. The frac- 
tion of the kinetic energy of the shells which is dissipated in the 
collision is 



/dyn — 



K+1 



(2) 



Then, if the injected kinetic power during the relativistic ejection 
phase is E, the Lorentz factor, comoving mass density and co- 
moving specific internal energy density in the shocked material 
are given by 

F, - — ^F, 

l+K 

E 



e* - 



47r/;2Fjc3 
2 



(3) 



These simple scaling laws will be used to explore the parameter 
space of the internal shock model in the next section. 

Once the dynamics of the internal shock phase is well un- 
derstood and the physical conditions in the shocked material are 
known, more assumptions are necessary to compute the emis- 
sion. This is described in the next subsection. 



2.2. Physical conditions in tine shocked medium 

The physics of the acceleration of particles in relativistic shock 
waves, as well as the amplification of the magnetic field, is far 
from being fully understood. It is therefore impossible in our 
state of knowledge to directly estimate the electron distribution 
and the magnetic field in the shocked medium from F,, and 
e, using first principles. Therefore, the microphysics related to 
these processes is usually parameterized in a very simple way, 
which is adopted in the present paper: (i) it is assumed that a 
fraction ee of the dissipated energy is injected in a fraction f 
of the ambient electrons that are accelerated to relativistic ener- 
gies, with a power-law distribution of slope -p. Note that most 
GRB studies (prompt and afterglow emission modelling) are re- 
stricted to the case ( - I (all electrons are accelerated) but nu- 
merical simulations of particle acceleration in relativistic shocks 
suggest that it may not be the case (see e.g. Byk ov & Meszaros 
1 1 9961: lEichler & Waxmanll2005t ISpitkovskvllloOSb ; (ii) it is as- 
sumed that a fraction es of the dissipated energy is injected in the 



magnetic field. We do not investigate in this paper an alternative 
scenario, where the magnetic field is dominated by a larg e-scale 
component anchored in the central source (see e.g. Spruit et aTl 
2001). With these four additional parameters (e^, p and e^), 
the number density of non-thermal electrons can be computed 



as well as their initial distribution 
n{y)^{p-\)-^\^\ forr>F„ 



with 



p - 2 fe P*e* p - 2 e^m^ e. 



(4) 



(5) 



(6) 



The magnetic field in the comoving frame of the shocked mate- 
rial is given by 



B' - VS^fBP.e* . (7) 
The evolution of F,n and B' is plotted for our example in Fig.[I] 

In practice, it is assumed that the relativistic electron dis- 
tribution extends up to a maximum Lorentz factor Fm, defined 
as the Lorentz factor where the acceleration timescale becomes 
comparable to the minimum of the radiative timescale and the 
escape timescale (see below). This corresponds to the most effi- 
cient acceleration that can be expected. In the comoving frame 
of the shocked region, the acceleration timescale of an electron 
with Lorentz factor y is estimated as R'^(y)lc, where 



<(r) = 



yntf-c 
eB' 



is the Larmor radius. This leads to 



Fm = min 



6ne 
a-jB' 



1/2 



eB'C 



(8) 



(9) 



where the radiative timescale is taken to be equal to the syn- 
chrotron timescale (Eq. ( fTTj i below) and the escape timescale 
is identified with the dynamical timescale f^x (Eq. ( flOl l below). 
Note that when inverse Compton losses are important, this ex- 
pression overestimates the maximum Lorentz factor Fm. This is 
further discussed later on. 



2.3. Emission in the comoving frame 

Timescales. Two timescales are necessary to characterize the 
physics in the shocked region: (i) the dynamical timescale 



R 

I> 



(10) 



which is the typical timescale associated with the adiabatic 
cooling due to the spherical expansion; and (ii) the radiative 
timescale f^^, defined as the timescale necessary for the rela- 
tivistic electron s to radiate most of their energy. As described in 
ISari et al electrons with t[.^^ <k f^^ are in "fast cooling" 

regime and will radiate efficiently, whereas electrons with f'^^ » 
are in "slow cooling" regime and will loose most of their 
energy via the adiabatic cooling. In internal shocks, the short 
variability timescale observed in the l ightcurves imposes that 
all electrons are in fast cooling regime jRees & MeszarosllI994l; 
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ISari et al.ll996HKobavashi et al.ll997h . This is also probably re- 
quired by pure energetic considerations, as the huge gamma-ray 
luminosities observed in GRBs are very difficult to understand 
if electrons are not radiating efficiently. From a numerical point 
of view, the advantage of being in fast cooling regime is that 
the emission is produced over a short timescale: relativistic elec- 
trons accelerated in one collision will radiate most of their en- 
ergy before the next collision occurs. This allows to compute the 
emission in an independent way: for each dynamical timestep 
(duration ~ f^^), the radiation in the shocked region is computed 
assuming that the dynamical quantities (e.g. the density) do not 
vary. 

Geometry. The shocked region is a shell with radius R, opening 
angle A6 (equal to the opening angle of the outflow, that can 
be considered as constant in the internal shock phase, the 
lateral expansion becoming efficient only when the outflow has 
notably decelerated) and comoving width A' ~ ct'^^ . During the 
dynamical timescale f^^, the emission in the comoving frame 
can be computed assuming constant dynamical quantities. The 
causally connected region during this duration has a size ct'^^ 
which is small compared to the lateral size of the shell RAO, 
as long as A9 » 1/f. In the comoving frame of the shocked 
region, one can therefore neglect the curvature and consider a 
infinite plane layer with width ct'^^. 

As electrons are in fast cooling regime with <K f^^, most 
of the evolution occurs on a short timescale, corresponding to 
a causally connected region of size cf^^ much smaller than the 
physical width of the region. Therefore, it is justified to assume 
that, if the shell is initially homogeneous, it will remain so for 
most of the evolution: the density distribution of electrons will 
depend on time, but not on the position in the shocked region. 
The same will happen for the photon distribution, which will 
appear as isotropic everywhere in the shocked region. This is of 
course not strictly valid within a distance cf,^^ from the edge of 
the shell, but the corresponding volume is negligible, as f^^ <K 



The photon field At time f - 0, just after the colHsion, when 
the particle acceleration and the amplification of the magnetic 
field are achieved (it is assumed that these processes operate 
on timescales which are short compared to the radiative and 
the dynamical timescales), the electron distribution is given 
by Eq. (|5]l and the photon density is zero. This is justified as 
all electrons that were shock-accelerated earlier have already 
cooled. 



an isotropic emission by electrons, this leads to 




= £ j| di'PAn, (13) 

where the homogeneity of the shock region is taken into account 
(see above). Finally, the photon density distribution is given by 

X'' P '(?) 

where it appears clearly that the local photon field is built by 
accumulating photons coming from a growing region of size ct' 
and therefore depends on the whole history of the emission be- 
tween ? = and P - t' . In the next paragraph this equation 
is expanded by giving explicitly the emission processes that are 
considered in the present study and including the absorption pro- 
cesses that were neglected in this paragraph. 

Radiative processes. Many radiative processes can operate in 
the shocked medium. In this paper, we focus on the processes 
that are expected to be dominant if the radiation is mainly 
produced by electrons, i.e. we do not include contributions 
associated to a possible population of relativistic protons 
accelerated in the sh ock. Such a component is included in the 
calculations made bv lAsano & In oue (2007) for a typical shock. 
Their results show that (i) for most parameters, the proton 
contribution is negligible, especially below a few GeV; (ii) it 
is only when » e^, i.e. when most of the dissipated energy 
is injected in the magnetic field and in protons, rather than in 
electrons, that a non negligible proton component emerges. 

Accelerated relativistic electrons in the amplified magnetic 
field will radiate via the synchrotron process. These synchrotron 
photons can be scattered to higher energies by relativistic 
electrons (inverse Compton). At low energy, they can also 
be absorbed (synchrotron-self absorption). At high energies, 
photon-photon annihilation can occur, producing electron- 
positron pairs. The corresponding pairs could contribute to the 
radiation, but this contribution is not considered in the present 
paper, as we limit our studies to cases where the production 
of pairs is weak (see next section). We di d not consider in 
this study the case of the "jitt er radiation" (iMedvedevI 120001 : 
Medvedev & Spitkovskvl |2008|) that is an alternative to the 
standard synchrotron radiation. 



The photon density distribution at time f at a given position 
in the shocked region is given by 



nv'(f') = 



c hv' 



(11) 



due to the isotropy of the photon field (see above). If absorption 
is neglected at this stage of the discussion, the specific intensity 
Iv'if ) is built by integration of 



dlyl 

ds' 



(12) 



Based on the timescales and the geometry discussed above, 
we have implemented a radiative code to solve the evolution 
of electrons and photons in the comoving frame of the shocked 
medium during a dynamical timestep. Two equations are solved, 
one for the evolution of the comoving electron density distribu- 
tion n (y, t') : 



dn 



d_ 
dy 



dy 
d? 



syn+ic+ad 



n (y, f) 



(15) 



and one for the evolution of the comoving photon density distri- 
bution riy'it') : 



along a ray, from i' = to s' - ct' (where i' = is the position 
where ny> is computed), due to the finite speed of light. Assuming 
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-CHy'it') Jdyn{y, t') cr^^{y,v') 
—cnyi(t') J" dv nyi(t')o-yy{v' ,v) 



(16) 



The indexes syn, ic, ad, sa and yy stand respectively for the 
following processes: synchrotron radiation, inverse Compton 
scattering, adiabatic cooling, synchrotron self-absorption and 
photon-photon annihilation. 

The expressions of the different terms appearing in Eqs. ( fTSt 
and ( fTSI l are listed in appendix |A] and the numerical method 
to solve this set of equations is described in appendix IB] 
The adiabatic losses are estimated by dy/dt'\^^ = ~T/'ex- 
The synchrotron radiation is computed exactly, assuming 
an isotropic distribution of the pitch angle a between the 
electron velocity and the magnetic field. The synchrotron self- 
abs orption is also compu t ed usi ng the exact cross-section (see 
e.g. iRvbicki & Lightrnanl Il979l) . Note that the corresponding 
heating term at low energy is neglected in Eq. ( fTSl ). Inverse 
Compton sc atterings ar e computed using the approximate kernel 
derived bv IJonesI (1968), which is a very good approximation, 
even in the Klein-Nishina regime. Note that Eq. ( fT6l ) does not 
include the loss term at low frequency corresponding to the 
source term at high energy. This is because the Thomson optical 
depth is always low in our case (see next section). We do not 
examine situations where comptonization could occur. Finally, 
the full cross-section for gamma-ga mma annihilation is used , 
assuming an isotropic photon field dGould & Schrede3ll967h . 
As mentioned above, the present version of the code does not 
include the associated pair creation term, so that we limit the 
study to cases where it is negligible (see next section). 



Following ISari et a n (ll998h . it is convenient to define Fn as 
the Lorentz factor of electrons whose synchrotron timescale 



^syn 



(r) 



6mneC 



is equal to the adiabatic cooling timescale, i.e. 



6nmeC 
O-T B' - fex 



(17) 



(18) 



When the synchrotron process is dominant, electrons with 
7 > Fc are in fast cooling regime. When inverse Compton 
scatterings become efficient, the effective transition between 
slow and fast cooling occurs at a Lorentz factor lower than F^ 
as the radiative timescale becomes shorter than the synchrotron 
timescale. 

The solution at t' - t'^^ of system of Eqs. dlST l and (fTST i is en- 
tirely determined by the expansion timescale f^^, the shape of the 
initial electron distribution (i.e. mainly Fm and p), the relativis- 
tic electron density n^'^'^ and the magnetic field B'. Rather than 
using these two last quantities, it is convenient to consider al- 
ternatively the critical Lorentz factor Fc and the initial Thomson 
optical depth associated to relativistic electrons 



4- = (Tjnl^'ct' 



(19) 



The radiative calculation has to be made for each collision 
occurring in the dynamical phase, i.e. at each instant along 
the propagation of a shock wave within the relativistic out- 
flow. Fig. 12] shows one of these elementary calculations. This 
case has been selected as the effect of each process is clearly 



identified (see caption of the figure). Possible additional ef- 
fects (scatterings or absorption) between photons emitted in 
a shocked region and electrons or photons present in another 
shocked region, which could affect the high-energy spectrum 
(iGruzinov & MeszaroslbOOO ). are not considered in the present 
paper but will be investigated in the future. We also ignore the 
effects of triplet pair production which can occur when electrons 
of very high energies encounter soft photons : the cross-section 
for this process becomes larger than the inverse Compton cross- 
section in the deep Klein-Nishina regime, for ye S 250, where 
y is the electron Lorentz fa ctor and e - hv/nieC^ the soft photon 
energy (Mastichiadis 1991). Finally, we do not include the possi- 
ble interaction of the prompt gamma-rays emitted from internal 
shocks and the circumburst environment, th at could also lead to 
an additional early high-energy component ( lBeloborodovl2002l 
2005.) . 



2.4. Observed flux 

Once the emission in the comoving frame is computed at each 
instant along the propagation of internal shocks within the rela- 
tivistic outflow, the observed flux as a function of time is com- 
puted by summing up the contributions of all shock waves, 
taking into account: (i) the relativistic effects (Lorentz trans- 
formation from the comoving frame of the shocked region to 
a fixed frame); (ii) the curvature of the emitting surface; (iii) 
the cosmological effects due to the redshift of the GRB source. 
The two first points require an integration over equal-arrival 
times surfaces, that is carried out following equations given in 
IWoods & Loebl (Il999l) . Any absorption in the gamma-ray range 
due to pair creation on the extragalactic background light is ne- 
glected. This would be important, depending on the redshift, 
above ~ lOGeV. Examples of synthetic lightcurves and spec- 
tra produced following the complete procedure described in this 
section are presented in Sect.|5] 

3. The emitted spectrum in the comoving frame 

As described in Sect.|2] the emitted spectrum in the comoving 
frame of the shocked material is entirely determined by four 
parameters: (i) the magnetic field B'; (ii) the adiabatic cooling 
timescale f^^; (iii) the relativistic electron density nl^'^, and 
(iv) the shape of the initial distribution of the Lorentz factor 
of accelerated electrons, i.e. the slope p and the minimum 
Lorentz factor Fm for a power-law distribution. A clear insight 
in the way that every of these parameters affects the radiative 
processes is necessary to anticipate the characteristics of the 
photon spectrum resulting from the shock-accelerated electrons. 

The final observed photon spectrum comprises the contribu- 
tions of all the photons emerging from the collisions occurring 
during the evolution of the relativistic outflow. We focus first on 
the radiative processes and the photon spectrum occurring after 
a single collision only and will describe later (§|5]l the complete 
GRB lightcurve and spectrum. 

We have carried out spectral calculations corresponding to 
a large exploration of the parameter space (B', t'^^, rif^, Fn,) de- 
scribing the physical conditions in the shocked material, as- 
suming a fixed electron slope p - 2.5. We computed 2744 
spectra corresponding to : (i) 7 values of the magnetic field 
log(B71 G) = 1, 1.5, 2, 2.5, 3, 3.5 and 4; (ii) 7 values of the 
dynamical timescale log(fg^/l s) = 0, 0.5, 1, 1.5, 2, 2.5 and 3; 
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Fig. 2. Emission in the comoving frame: an example. This figure shows the result of the radiative calculation in the comoving 
frame of the shocked material, at time 1.9x10^* s in the example shown in Fig. [1] Left. Evolution of the relativistic electron population. 
The initial distribution is plotted with a thick line. Right. Evolution of the photon spectrum v'^Wy'- The final spectrum at t' - t'^^ 
is plotted with a thick line. Several phases are well identified: the synchrotron spectrum is built very early (thin solid line). Once 
these seed photons are present, inverse Compton scatterings produce high energy emission (dotted line). At late times, high-energy 
photons annihilate which reduces the intensity of the second peak in the spectrum (dashed line). The microphysics parameters are 
fe = 1 /3, ^ = 10 p -2.5 and es = 1 /3 so that the initial Thomson optical depth is Tj'^ ^ 6 x 10"^ and the critical Lorentz factor 
equals Fc - 2.4 (see text). 



(iii) 7 values of the electron density log («e"/l cm ^) = 4, 5, 6, 
7, 8, 9 and 10; (iv) and 8 values of the minimum electron Lorentz 
factor log F,n = 1 .5, 2, 2.5, 3, 3.5, 4, 4.5 and 5. 

3.1. Radiative efficiency and transparency 

In the shocked medium, the evolution of the relativistic electron 
distribution is governed by several radiative processes that are 
in competition with the adiabatic cooling due to the spherical 
expansion. The efficiency of converting the energy deposited in 
relativistic electrons in radiation depends strongly on the rela- 
tive magnitudes of the radiative cooling timescale t'^^^ of rela- 
tivistic electrons and the adiabatic cooling timescale of expand- 
ing shell (see Sect. |2.3l l. The observed short timescale vari- 
ability as well as the high isotropic equivalent energy radiated 
in gamma-rays imply that electrons are radiati ng efficiently, i.e . 
that f^^^ < fg^, the so-called fast-cooling regime (ISarietal.119981) . 
Therefore, in the following, we have only considered the region 
of the parameter space where the radiative efficiency is high, i.e. 
Umd/uf'^ > 0.5, where 

Mg"^*^ — J" dy n (y, t' = 0) ym^c^ 

is the initial energy density in relativistic electrons and 

Miad - ^ dv' riy' (t' = fg^) hv' 

is the final energy density contained in the radiated photons. 
Fig.|3]shows this region in the plane Fm/Fc-Mic/Msyn, where Msyn 



(resp. Mic) is the component of Miad corresponding to synchrotron 
emission (resp. inverse Compton emission). Clearly, when 
inverse Compton scatterings are inefficient, the electron radia- 
tive timescale is the synchrotron ti mescale and the efficiency 
condition is equivalent to F^ > Fc dSari et al.lll998h . However, 
inverse Compton scatterings reduce the effective electron radia- 
tive timescale and models with F^ < Fg can still be efficient if 
Mic > Msyn. In these models, the Lorentz factor fg of electrons 
whose radiative timescale is equal to t'^^ can be much lower than 
Fc (defined from the synchrotron timescale only). 

As we do not consider scenarios where the emitting re- 
gion is opti cally thick (for instance a comptonized spec- 
trum, see e.g. Ghisellini & Celottill999HMeszaros & Reesl2000l; 
Pe'er & Waxman 2004), we also limit the discussion to the re- 
gion of the parameter space where the medium is optically thin 
for Thomson scatterings, i.e. tt - o-jn^ct'^^ < 0.1, being 
the total density of electrons (relativistic or not). Using the two- 
shells model, this condition leads to a minimum value for the 
Lorentz factor: 



F > 1.58 



64/^7/2 



Anm^c'^T 



1/5 



(20) 



This minimum Lorentz factor of the outflow is of the order of 
100 - 200, for T ^ Is, £■ ^ lO'^^ erg s"' and ^ 2 - 4. However, 
an additional effect must be taken into account: high energy pho- 
tons can annihilate into e"^ pairs, and the corresponding new lep- 
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Fig. 3. Radiative efficiency. The explored region of the param- 
eter space of internal shocks is shown in the plane Mic/Msyn 
versus Fm/rc (shaded region labeled 'All'). The x-axis quan- 
tity measures the efficiency of inverse Compton scatterings and 
the y-axis compares the synchrotron and the adiabatic cooling 
timescales. Radiatively efficient cases, i.e. cases where most 
of the energy initially injected in relativistic electrons is radi- 
ated (see text) are indicated (shaded region labeled 'Rad. effi- 
cient'). Inverse Compton scatterings reduce the effective radia- 
tive timescale of relativistic electrons. Therefore models with 
large Mic/Msyn Can be efficient even for low Fm/Tc. For negligi- 
ble inverse Compton emission, the radiatively efficient region 
corresponds to the usual synchrotron fast cooling regime. 



tons will increase the Thomson optical depth. Therefore the true 
transparency condition that we impose is 



T!r°' = TT(l+^«,/«r)<0.1, 



(21) 



where n± is the final density of leptons produced by pair 
annihilation. This will increase the minimum value of the 
Lorentz factor derived above. Compared to an alytical estimates 
of th e minimum Lorentz factor (see e.g. Li th wick & Sm| 
1200 lb . we use here a precise estimate of the pair production 
factor n±lnf''^ which is a byproduct of our radiative calcu- 
lation. ' Lithwick & Saril (l2001i) have also considered a third 
transparency condition to derive a minimum value for the 
Lorentz factor from GRB observations: the absence in the MeV 
spectrum of a cutoff due to yy annihilation. We will discuss this 
condition in Sect. 14.41 as it is expected that this cutoff could be 
observed by Fermi in some GRBs in the future. A consequence 
of our transparency condition is that all the cases presented in 
this paper correspond to situations where the fraction of energy 
(initially in high energy photons) which is deposited in pairs 
is small. This justifies that the emission of these leptons is 
neglected in our present calculations. 

Except for the two limitations (efficiency and transparency), 
all parameters are a priori acceptable. Indeed, models for the 



central engine of gamma-ray bursts are not in a state where a 
distribution function can be provided for the injected kinetic 
power or the initial Lorentz factor in the outflow. Even the 
expected range of each quantity is highly uncertain. Therefore, 
the most promising way to estimate such physical quantities is to 
apply a detailed spectral model as described in this paper to re- 
cover the internal shock parameters that can reproduce observed 
lightcurves and spectra. When the low energy gamma-ray range 
only (e.g. BATSE data) is used, there is a large degeneracy. 
Hopefully, observations in the high-energy gamma-ray range 
{Fermi data) will improve this situation. For this reason, we 
study in this section how the broad spectral shape is affected by 
each parameter of the model. Before this, we recall the main 
scaling laws that are expected from analytical considerations, 
and check their validity with our detailed calculation. 



3.2. Analytical estimates 

Synchrotron component. The dimensionless photon frequency 
in the comoving frame is defined by e = hv' jm^c^. The exact 
solution of Eqs. (fTsT i and (fTST i can be obtained when only syn- 
chrotron radiation and adiabatic cooling are included. The time- 
averaged (over fg^) electron distribution h{y) is very close to a 
broken power-law: 

1 I (f)"^ ffFc <r<rm. 



n(r) 



(i^r(x)-"'-' if,>i, 



(Fn, < Fc) . 



(22) 

An accurate approximation of the corresponding ph oton spec- 
trum Wflsyn is given by a broken powerlaw shape jSari et al.l 
Il998h : 



~ (tr^ if.c<f<.™, (23) 

(^f if e > em , 
in the synchrotron fast cooling regime (F^ > F^) and 

(^)^ ^ ifem<e<ee, (24) 

in the synchrotron slow cooling regime (F„i < Fc). In these 
expressions, (resp. e^) is the synchrotron frequency of elec- 
trons at Fm (resp. Fc). It appears clearly that the first case (syn- 
chrotron fast cooling regime) is efficient as most of the initial 
energy deposited in relativistic electrons is radiated. In the sec- 
ond case (synchrotron slow cooling regime), only a small frac- 
tion of the electron energy is radiated. The peak of the syn- 
chrotron emission in e^n^ (equivalent to vFy) is located at fre- 
quency£p_syn = max {en,; ed for 2 < < 3 (em if > 3), i.e. from 



\4/3 



Isyn 



P-2 



Eq. ( IA.4l i at energy 



hy'^,^. - 0.17 eV 



B' 



1000 G 



^ p 
100 



with Fn 



max {F„ 

Fm 



if 2 < /:>< 3 
if/? > 3 



(25) 
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Our estimate of the maximum Lorentz factor of relativistic elec- 
trons (Eq. (|9]l) leads to a cutoff in the synchrotron component at 
frequency 



/zv;,.,,„-min^230MeV; 5300 TeV j ^ ^) ' 



(26) 



Except for very low magnetic fields and very short dynamical 
timescales, it is always the first limit, when the acceleration 
timescale becomes of the order of the synchrotron timescale, 
that dominates. 

As shown in Fig.|4] (bottom panels) and Fig.|5] (case (a)), 
the numerical results of our radiative code when synchrotron 
radiation is dominant are in an excellent agreement with 
these analytical estimates, showin g the goo d accu racy of the 
synchrotron spectrum described bv lSari et al] (Il998l) . 



The timescale associated with the synchrotron self- 
absorption at frequency v' can be estimated by 



2- 



d\nn 



lilny 



(27) 



From Eq. (l22l i. one gets in the synchrotron fast cooling regime 



t'Ay') 47rv- 



if v' < , 
if v' < v' < v' 



(28) 



, (yjy'c) (.y'/y'm) - '^^y' > y'm 



and in the synchrotron slow cooling regime 

;;+4 



f'(v') Anv''V, 



if v' < v' < v' 



(29) 



At high frequency, this timescale is very long and the syn- 
chrotron self-absorption process is negligible. It will only affect 
the spectrum below frequency for which t'jy.^ - t'^^. One can 
show that below v^, the predicted slope of the absorbed spec- 
trum is e^Hj oc €^ if < v^, and e^n^ oc e^^^ otherwise. We 
checked with our numerical code that the accuracy of these ex- 
pressions is quite good as long as the inverse Compton cooling 
is not dominant (othe rwise the true e lectron distribution differs 
from Eq. ^ used bv lSari et al] (Il998l) . see below). 

Inverse Compton component. If most scatterings between rel- 
ativistic electrons and synchrotron photons occur in Thomson 
regime, the peak of the inverse Compton component is expected 
^ic,p ~ -Fp6p syi-i, i.e. 



W,^^ 1.7keV (— ^UIl 
P'"= 1 1000 G/ 1 100 



(30) 



The Thomson approximation is valid as long as Fpfp syn ^ 1, i e. 



Tp « Fkn ^ 3.0 X 10' 
which corresponds to 



B' 



1000 G 



-1/3 



^y'^,ic ^ ^^KN - rKN»ieC^ - 1-3 GeV 



B' 



1000 G 



'1/3 



(31) 



(32) 



A severe reduction of the high-energy spectrum should always 
be expected above vj^j^ in the comoving frame due to Klein- 
Nishina corrections to the inverse Compton cross-section. Even 
more, the maximum Lorentz factor of the relativistic electrons 
given by Eq. (|9]l leads to an absolute maximum energy ~ FmOTcC^ 
for a scattered photon, i.e. 



B' 

hvi. ■ ^ min <! 1 -9 Te V ! 

M,ic 1 'lOOOG 



-1/2 

; 9000 TeV 



B' 



1000G/\1 s 



(33) 



Again, except for very weak magnetic fields and very short dy- 
namical timescales, the maximum inverse Compton frequency 
is always given by the first limit (acceleration limitation due to 
radiative losses). From these estimates, one can deduce that the 
peak of the inverse Compton component should be found in all 
cases at the frequency v' ^ minjvpj^; v^^^; vj^j^j. Fig. |4] (top 
panels) show that numerical results of our radiative code are in a 
reasonable agreement with these theoretical predictions, as long 
as inverse Compton scatterings are not the dominant cooling 
process for electrons. The scaling in equation ( l30l l appears to 
be correct but the normalization seems to be underestimated 
by a factor ~ 10. When Klein-Nishina corrections become 
important, the peak of the inverse Compton component appears 
typically at frequency ~ O.lvJ^j^ compared to Eq. (|32]|. This is 
not surprising since the cross section for scatterings of a photon 
with frequency v' by an electron with Lorentz factor y shows 
non negligible Klein-Nishina deviations well below the limit 
hv' = ym^c^. Significant deviations of one order of magnitude 
from the simple estimate 0.1 vj^ can be observed (see the cases 
with a weak magnetic field in Fig.|4|l. 

When the Thomson regime is valid, the ratio of the inverse 
Compton over the synchrotron power is given by the Compton 
parameter, defined as the ratio of the energy density in pho- 
tons over the magnetic energy density, Y - Mnd/ (b'^ /Stt). 
This quantity is time-dependent. However, when not stated oth- 
erwise, Y stands in this paper for the final value of the Compton 
parameter at f^x- As long as F < 1, synchrotron losses dom- 
inate, the seed photons for inverse Compton scatterings have 
the spectrum which is given above by Eqs. ( |23] l and (l24l i. and 
the distribution of the electrons responsible for the scatterings 
is close to the brok en-power law distribution n(y) described 
by ISari et all (Il998l) . The corresponding spectral s hape of the 
invers e Compton component has been derived by Sari & Esiiil 
(I2001h and is given in their appendix A. It is based on the inte- 
gration of the approximate relation 



dy n(y) 



hv' 



(34) 



where h(y) is the time-averaged electron distribution predicted 
by Eq. (|22] | and PyJiy) is the inverse Compton power radiated 
at frequency v' by an electron with Lorentz factor y, computed 
assuming a seed photon distribution equal to the standard, time- 
averaged, synchrotron spectrum given by Sari et al. (1998]). 
Instead of the complete expression of P'^Xy) (see appendix [All, 
the authors use a simplified kernel (which is equivalent to 
assume Thomson regime everywhere) so that the integration 
can be made analytically. In the present study, both the electron 
and photon distributions are time-dependant, which leads to 
significant differences in the high-energy spectrum compared to 
the time-averaged approach. This will be further discussed later 
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Fig. 4. Synchrotron and inverse Compton peaks: comparison with analytical estimates. The photon peak energy of the syn- 
chrotron (bottom) and inverse Compton (top) components is plotted as a function of the initial minimum Lorentz factor of the 
accelerated electrons Tm, for five different values of the magnetic field (color coded) and increasing values of the accelerated elec- 
tron density n^'^'^ in the three panels. The dynamical timescale is fixed to t'^^ = 10 s. Values obtained from our radiative code are 
indicated with open (resp. filled) circles corresponding to calculations without (resp. with) yy annihilation. In the bottom panels, 
when the circle is filled with a black cross, the inverse Compton process is dominant compared to the synchrotron radiation. These 
results are compared to the analytical estimates (solid lines) discussed in the text: Eq. dZSl l for the synchrotron component and 
Eqs. ( [30l ) renormalized by a factor 10, ( [32l ) renormalized by a factor 0.1 and (33[ for the inverse Compton component. Note that, in 
the bottom panels, the transition between the horizontal and rising branches of the analytical estimate for fp syn corresponds to the 
transition between the synchrotron slow and fast cooling regimes. 



on. 



The intensity of the inverse Compton component is 



Isyn 



(35) 



The Compton parameter in this case equal^ 



1 



3p-2 



1 ^ ifr„, >r,, 

{y^Y ' if < fc and 2 < p < 3 , 
Tm/Fc if Fn, < Fc and p >3 . 

(36) 



Note that the term ^j^r^'^T^r^ is simply equal to u\ 



/mb, i.e. 

fe/fB using the standard parameterization of the microphysics. 
In synchrotron fast cooling regime, the value of Y can be un- 
derstood as follows: due to a short synchrotron timescale, the 
region populated with relativistic electrons has a typical size 
~ cfjy^ (Fm). The effective Thomson optical depth for relativis- 
tic electrons is ~ o-jnl^'^ct'^y^^ (F,n). In Thomson regime, the fre- 
quency of a photon scattered by an electron with Lorentz factor 
y is multiplied by y^. Therefore, the Compton parameter Y is 
approximatively given by 

Y-^^hnT'ct:yJTj)ri (37) 

' These expressions assume that the maximum electron Lorentz fac- 
tor Fm is greater than Fc, which is always the case in the fast cooling 
regime. On the other hand, a "very slow" cooling regime is possible 
when F,^ < Fm < F^. In this case the break at in the synchrotron 
spectrum is suppressed as it is above the cutoff at v^syn- Compton 

parameter in this case equals ^ - ^ {t^) (f^) ' if 2 < p < 3 and 



y.|(r^)ifp>3. 



which leads to the expression given above. On the other hand, 
in synchrotron slow cooling regime, the size of the region 
populated by relativistic electrons is now given by ~ c^x and 
therefore the Thomson optical depth for relativistic electron 
becomes much larger ~ cr-Ynl^'^ct'^^, which leads to Eq. (|36l ) for 
p > 3. In the case where 2 < p < 3, this equation is coiTected 
by a factor (Fc/Fm)^"'' > 1 as the typical Lorentz factor of 
scattering electrons is at an intermediate value between Fm and 
Tc. 

As shown in case (a) of Fig.|5] except at very high energies 
(above /ivj^j^), our numerical calculations are in very good 
a greem ent with the approximate spectrum given bv lSari & EsinI 



200 ih . as long as F < 1. Above v'^^, a better estimate 
obtained when integr ating numeri c ally E g. ( l34l i with the same 
assumptions used by ISari & EsTnl (1200 ll) (time-averaged elec- 
tron and seed photon distributions) but using a complete kernel 
that includes Klein-Nishina corrections in the inverse Compton 
power (thin black line in Fig.|5]). However, even using a more 
accurate cross section, Eq. ( [34l ) always overpredicts the inverse 
Compton emission at high energy. This systematic difference 
appears because the high energy photons in the inverse Compton 
component are due to the scatterings of photons at v' ~ vj, 
with high Lorentz factor electrons at y > max{Fin;Fc). In fast 
cooling regime, these two species are not present at the same 
time in the shocked region, as the duration necessary to form the 
synchrotron spectrum at Vp is also the duration necessary to 
cool electrons above Fm, i.e. the synchrotron timescale f^^^ (Fm). 
As Eq. ( [34l i is based on a time-averaged approach, it cannot 
take into account such effects, related to the way the radiation 
field is built, and therefore it overestimates the spectrum above 
Vpj^. We checked that in synchrotron slow cooling regime the 
agreement is better above v^j^ than what is observed in Fig.|5] 




Fig. 5. Synchrotron and Inverse Compton spectral components: comparison with analytical estimates. The time-averaged 
electron distribution and the final photon spectrum are shown for 3 different cases: (a) r,,, - 10^ and n^'^'^ - 10^ cm"^; (b) = 10^ 
and nf^ = 10'° cm"^; (c) = lO"* and nf^ = 10** cm^l In all cases, p = 2.5, B' = 3000 G and 4 = 10 s. Synchroti-on 
self-absorption and yy annihilation are not included in the calculation. The time-averaged electron distribution obtained with our 
radiative code is plotted in the top panel with a red line and compared to the analytical prediction in black. The characteristic 
electron Lorentz factors discussed in the text are indicated with vertical bars. The initial distribution is also plotted with a dotted 
red line. The synchrotron (resp. inverse Compton) component of the final photon spectrum is plotted i n the bottom panel with a red 
(resp. blue) line and is compared to the analytical estimate in black (thick line: analytical estimate bv lSari & EsinI |2001); thin line: 
semi-analytical estimate described in the text, which includes Klein-Nishina corrections). For indication, the result of the numerical 
calculation including yy annihilation is also shown in dashed red and blue lines. The characteristic photon frequencies discussed 
in the text are indicated with vertical bars. Klein-Nishina corrections are negligible (except at high energy) in case (a) and (b). 
Synchrotron radiation is dominant in cases (a) and (c), whereas inverse Compton scatterings dominate in case (b). 



but Eq. ( [34] l is still overestimating the high-energy component 
when evaluating the scatterings by fast cooling electrons (i.e. 
electrons with y > Fc). 

When F > 1, inverse Compton losses become dominant. 
Then, the effective radiative timescale is shorter than the syn- 
chrotron timescale (by a factor ~ (1 -H Y)), the effective critical 
Lorentz factor is reduced (f ^ ~ F^ /( 1 -H F)) and the corresponding 
frequency ec in the synchrotron spectrum (Eqs. ( l23l l and ( l24b ) as 
well. The intensity of the synchrotron component is reduced by a 
f actor 1 /( 1 + y) . These corrections are however very approximate 
and valid only in Thomson regime. Our tests show that when in- 
verse Compton scatterings become dominant, the modified cool- 
ing rate of electrons affects the time-averaged distribution n(y) 
(which differs from the standard broken-law distribution given 
by Eq. (|22])). and therefore the distribution of seed synchrotron 
photons becomes different from the standard synchrotron spec- 
trum given by Eqs. ( l23T l and ( |24] |. This is well seen in cases 
(b) and (c) i n Fig. [31 In fact, in this case, the approach used by 
ISari & EsinI (|20()1 ) is not appropriate because the spectrum of 
the seed photons cannot be predicted by an a priori calculation 
including the synchrotron process only: the resulting spectrum 
has not enough time to be built when inverse Compton losses 
are included. This effect becomes stronger when Klein-Nishina 
corrections are important, as the ratio of the inverse Compton to 
the synchrotron power becomes highly dependant on the elec- 
tron Lorentz factor. As seen in Fig.|5] the low-energy slope of 
the synchrotron spectrum is steeper in that case. Such a be- 
havior is in agreemen t with the theoretical predictions made by 
iDerishev et alj ( l200ll) . We plan to investigate in a forthcoming 
paper if this could reconcile the synchrotron radiation with the 



observ ed distribution of th e low-energy photon index in BATSE 
bursts ("P reece et al.ll200 C^. which differs from the simpl est pre- 
diction of the fast cooling synchrotron spectrum CGhisellini et al.l 
l2000h as its mean value is close to o- L 



Formation of the radiation field. These results show that the 
high energy component of the photon spectrum cannot be 
estimated accurately without understanding how the radiation 
field (seed photons for inverse Compton scatterings) is formed. 
Initially, no photons are present and synchrotron radiation is 
always dominant {Y(t' = 0) = 0). The Compton parameter is 
an increasing function of time, due to the progressive building 
of the radiation field (see Fig. lC.lb . When synchrotron radiation 
is the dominant process, the radiation field increases up to 
~ ^syn(rm) - 4x^0 /Tm in fast cooling regime, and then it 
saturates. In slow cooling regime, it increases up to f ^ t'^^. 
In both cases, the time evolution of the Compton parameter, 
y(f'), can be evaluated analytically (as shown in appendix ICl) 
and its asymptotic value Fex = Y (f^^) is given by Eq. (|36] |. 
as long as most scatterings occur in Thomson regime. In fast 
cooling regime, a necessary condition to have a dominant 
inverse Compton component in the final spectrum is therefore 
uf^ Iub - Se/eB > 1. It is however not a sufficient condition, 
as the intensity of the inverse Compton component can be 
attenuated by Klein-Nishina effects, and also by yy annihilation. 

When Fex > 1, the impact of inverse Compton scatterings 
on the electron distribution will depend on the time f'^ where 
y(f'^) ~ 1, i.e. the time when inverse Compton scatterings 
become the dominant process of cooling. Indeed, only the 
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distribution of electrons below rc(f[^) can be affected by the new 
dominant cooling process, as electrons at higher Lorentz factor 
have already cooled by synchrotron radiation. Here, the Lorentz 
factor Tcit') is defined as the Lorentz factor giving a synchrotron 
timescale of the order of t\ i.e. fsy„(rc(f')) - t'. With this 
definition Fc = rdt'^^). In the synchrotron fast cooling case, the 
synchrotron spectrum around the peak at will be affected by 
inverse Compton scatterings if this process becomes dominant 
at very early times, i.e if ^dtU > ^m, which is equivalent to 

When inverse Compton scatterings are extremely efficient, 
they can represent the dominant electron cooling process, even 
at early times. When this happens, the maximum Lorentz factor 
of accelerated electrons Fm is overestimated in Eq. (|9]l. From the 
evolution of Y{t') discussed in appendix|C] one can deduce the 
value reached by the Compton parameter when electrons at Fm 
have cooled, i.e. 



(38) 

Here it is assumed that t'^^^ (Fm) <k fex Eq. (IC.4l l. which is al- 
ways the case when the radiative cooling is the dominant limiting 
process for electron acceleration. If I'm < 1, the estimate of the 
electron maximum Lorentz factor Fm given by Eq. (|9]l is valid as 
the radiative timescale of electrons at Fm is accurately given by 
their synchrotron timescale. On the other hand, if Fm > 1, it is 
possible that the value of Fm given by Eq. ^ is overestimated. It 
is not always the case as the value of Ym is computed assuming 
inverse Compton scatterings in Thomson regime. The true value 
of Fm can therefore be reduced by Klein-Nishina corrections. In 
practice, we checked that our assumptions regarding the maxi- 
mum Lorentz factor of electrons are consistent in all the cases 
presented in this paper. 

Photon-photon annihilation. The timescale associated with yy 
annihilation is given by 



dv' 



5/! 



1 



ccTj {hv'riy') 



,(39) 



where a Dirac approxim ation has been used for the cross section 
jGould & Schrede3ll96 7). The cutoff will occur at high energy 
and the corresponding photons will annihilate with low-energy 
photons whose distribution is approximatively given by the syn- 
chrotron spectrum described in Eqs. ( 1231 ) and ( l24l i. An approxi- 
mate shape of the absorbed spectrum can then be computed by 
attenuating the emitted spectrum by a factor 



t' (■/) 

!zzllZ(i_e-':,/';,(v')) ^ 



(40) 



where Tyy(v') is the 77 annihilation optical depth at frequency 
v'. A comparison with the results of the detailed radiative code 
shows that this approximate treatment is again accurate as long 
as inverse Compton losses are not the dominant cooling process 
for electrons, i.e. as long as the low-energy photon distribution is 
well described by the standard synchrotron spectrum. In Fig.|5] 
the attenuation of the spectrum at high energy due to 77 annihi- 
lation is shown in different cases with an increasing importance 
of inverse Compton scatterings. 



3.3. The shape of the radiated spectrum 

We define a "reference case" corresponding to the physical 
conditions in the shocked material of the example in Fig.[T] at 
f = 1.9 X 10^ s, i.e. F, =^ 240, p. ^ 6.8 x 10"'^^ g cm^^ 
ejc^ ^ 8.0 X IQ-^ and 4 = 80 s. For = eg = 1/3, = 2.5 
and ^ = 10"^ (only 1 % of the electrons are accelerated), this 
leads to B' ^ 2000 G, nf" ^ 4.1 x 10^ cm"^ and F„ = 1.6 x lOl 
This choice of parameters is motivated by the study presented 
in Daigne & MochkovitchI ([1998.) which favors the case where 
the magnetic field is high and where only a small fraction of 
electrons is accelerated, as these two conditions are required for 
the synchrotron peak to be in the BATSE range. Starting from 
this "reference case", one of the parameters (Fm, B', t'^^, n'^'^) is 
varied, while all other parameters are maintained constant. The 
resulting evolution of the spectrum is plotted in Fig.|6] 



Effect of the initial minimum Lorentz factor of relativistic elec- 
trons. Panel (a) shows the effect of F,n. The spectrum is the 
combination of a low-energy component due to the synchrotron 
emission in fast cooling regime, and a high-energy component 
due to inverse Compton emission, partially suppressed due to 77 
annihilation, (i) Synchrotron emission. The evolution of the low- 
energy peak (synchrotron emission) with Fm follows exactly the 
prediction of the analytical estimate in fast cooling regime: the 
spectrum peaks at the frequency v'^ (synchrotron frequency of 
electrons with 7 - Fm) that scales as oc F^^ and the corre- 
sponding peak intensity follows {v'^nyi^ , Fm (ISari et al.l 

119981) . A dashed line of slope 1 /2 in the log v'-log ( v' ^ Wv') dia- 
gram indicates the predicted position of the synchrotron peak. 
The agreement with the numerical calculation is excellent; (ii) 
inverse Compton scattering. For low Fm, as synchrotron pho- 
tons peak at low energy, inverse Compton scatterings occur in 
Thomson regime. Then, from Eq. ( |36] |. the Compton parameter 
Y scales as F oc Fm. This is indeed shown in panel (a), where the 
intensity of the inverse Compton component increases when Fm 
increases. The high-energy peak due to inverse Compton emis- 
sion follows exactly the predicted line of slope 1 /2 in the log v'- 
log {v'^ Hy') diagram (the peak energy scales as Fm Vm oc Fm and 

the intensity scales as F^v'^ny') ^ Tm)- Fm increases, the 
synchrotron emission peaks at higher energy and more and more 
photons have energies comparable with m^c^ in the frame of 
electrons at 7 = Fm. Therefore the efficiency of inverse Compton 
emission is strongly reduced by Klein-Nishina corrections; (iii) 
77 annihilation, this process becomes more efficient for high Fm 
as the synchrotron emission peaks at higher energy and more 
photons are above the threshold. In conclusion, these two effects 
combine so that the high-energy component is most intense for 
intermediate values of Fm. 



Effect of the magnetic field. Panel (b) shows the effect of B' . 
This effect is more complicated than for Fm, especially at low 
energy. ( /) Synchrotron emission. For very weak magnetic fields, 
inverse Compton scatterings are so efficient that the effective 
electron radiative timescale is short compared to the adiabatic 
cooling time and most of their energy is radiated. This high ra- 
diative efficiency is reached despite the fact that synchrotron ra- 
diation occurs in slow cooling regime {t[^^ (Fm) > t'^^). In other 
terms, this situation corresponds to the following ordering of 
characteristic electron Lorentz factors: Fc < Fm < Fc. For more 
intense magnetic fields, the electron synchrotron timescale is re- 
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Fig. 6. Emitted spectrum in the comoving frame. We consider a "reference case" defined by - 1600, B' - 2000 G, 4^ = 80 s 
and n'f^'^ =4.1 x 10^ cm"^ corresponding to the physical conditions in the shocked material at f = 1.9 x lO'* s in the example shown 
in Fig.[T] for - - ^ /3, p = 2.5 and ^ = 10"^ (see text). Each panel shows the evolution of the spectrum when one parameter 
is varied, while all other parameters are maintained constant: (a) effect of the initial minimum Lorentz factor of the electrons Fn,. 
For low Fn,, the weak inverse Compton component is shown with a dotted line. The high-energy component that would be obtained 
without yy annihilation is plotted with a thin solid line. Two dashed lines indicate the predicted position of the synchrotron peak in 
fast cooling regime, and of the inverse Compton peak in Thomson regime; (b) effect of the magnetic field B'. For low B', the weak 
synchrotron component is plotted with a dotted line. The high-energy component that would be obtained without yy annihilation 
is plotted with a dotted line for B' = 6.3 G and B' = 2000 G. The dashed line indicates the predicted position of the synchrotron 
peak in slow and fast cooling regime. ; (c) effect of the adiabatic cooling timescale f^x- For t'^^ = 80 s, the spectrum that would be 
obtained without yy annihilation is plotted with a dotted line; (d) effect of the initial density of relativistic electrons nl'^^. For high 
nl'^'^, the spectrum that would be obtained without yy annihilation is plotted with a dotted line. In all panels, the approximate range 
of low-energy gamma-rays (1 keV-1 MeV) is indicated by the gray area, assuming z = 1 and F* = 240. 
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duced and synchrotron radiation operates in fast cooling regime 
(fc < Tc < Fn,). The complex evolution of the synchrotron 
peak with B' observed in panel (b) is then due to the transition 
between synchrotron slow (low B') and fast (high B') cooling 
regime. For low B', the synchrotron peak is given by the syn- 
chrotron frequency of electrons with y = Fc and decreases 
with B' as oc B' ^^. The corresponding peak intensity follows 
(v'^Mv')^, ^, oc b2(p-2) dSari et al.llT998i) . so that the low-energy 

peakfollows a line of slope -2(/7-2)/3 in the log v'-log ( v' ^ My') 
diagram. For larger values of B', synchrotron radiation is in fast 
cooling regime and the peak is given by v'^, which scales as 
v'^oz B'. The corresponding peak intensity does not vary with B' 
and therefore the low-energy peak falls into an horizontal line for 
high magnetic fields in the same diagram; (ii) inverse Compton 
scattering. For low magnetic fields, when the synchrotron radia- 
tion is in slow cooling regime, the typical size of the region pop- 
ulated by relativistic electrons is ~ ct'^^, that is much larger than 
in the synchrotron fast cooling regime. Therefore, even in Klein- 
Nishina regime, inverse Compton scatterings are extremely ef- 
ficient due to the increased optical depth. As B' increases, the 
synchrotron peak energy decreases and inverse Compton scatter- 
ings can occur in the Thomson regime. However, as soon as the 
magnetic field is high enough so that the synchrotron radiation 
enters the fast cooling regime, the Thomson optical depth for rel- 
ativistic electrons is reduced as the size of the region populated 
with relativistic electrons is ~ ct[^^ (Fm) <K cf^x- In addition, the 
synchrotron peak energy increases, so that Klein-Nishina correc- 
tions are present. These two effects result in a strong decrease of 
the inverse Compton efficiency; (iii) yy annihilation. This pro- 
cess becomes important for a low magnetic field B' , when the 
emission at high energy is most efficient. In conclusion, weak 
magnetic fields favor a strong high-energy emission, even when 
yy annihilation is non negligible. 



Effect of the adiabatic cooling timescale. The effect of t'^^ is 
shown in panel (c). For a magnetic field B' - 2000 G and an 
initial minimum electron Lorentz factor Fn, = 1600 as in our ref- 
erence case, the synchrotron timescale equals f^^^ (Fm) - 0.12 s. 
(i) for high values, f^x » 0.1 s, the spectrum at low energy does 
not depend on t'^^ as it is produced by synchrotron radiation in 
fast cooling regime. Inverse Compton scatterings are rare be- 
cause relativistic electrons are present only in a small fraction 
of the volume of the shocked region (~ fsyn/fex' see above). On 
the other hand the absorption due to yy annihilation increases 
with fgx because the effective size of the region populated with 
photons is larger and the probability to have photon-photon in- 
teractions is increased; (ii) for low values, t'^^ <K 0. 1 s, the radia- 
tive efficiency decreases strongly as electrons are in slow cooling 
regime; (iii) for intermediate values such as t'^^ ^ f^y^ (Fm), rela- 
tivistic electrons can radiate efficiently and still populate most of 
the shocked region. This favors a brighter high energy compo- 
nent as inverse Compton scatterings are more frequent. In con- 
clusion, an intense high energy component in the photon spec- 
trum is favored by intermediate values of the adiabatic cooling 
timescale. 



Effect of the initial density of relativistic electrons. The effect 
of rif^ is shown in panel (d). At low densities, the optical depth 
for inverse Compton scatterings is low, so the spectrum is simply 
a synchrotron spectrum in fast cooling regime. At high densities 
the intensity of the inverse Compton component increases. 
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Fig. 7. Physical conditions in the shocked region (comoving 
frame) that favor an intense high energy component. From 
our exploration of the parameter study of the internal shock 
model (see Sect.|3]l the histograms of the electron minimum 
Lorentz factor F^, the magnetic field B' , the adiabatic cooling 
timescale t'^^ and the initial density of accelerated electrons nf^ 
are plotted for all models that satisfy the radiative efficiency and 
transparency conditions (see text). The shaded histograms show 
the models that fulfill the same conditions and, in addition, have 
an efficient inverse Compton emission (most of 50% of the elec- 
tron energy which is radiated is due to inverse Compton scatter- 
ings) and a negligible yy annihilation (defined as cases where 
the peak in vFy of the high-energy component is shifted by less 
than 10 % when including yy annihilation). 



as well as the absorption due to yy annihilation, because the 
number of available high-energy photons scales with the density 
of emitting electrons. Therefore the most intense high-energy 
component is again obtained for intermediate values of the 
relativistic electron density, when inverse Compton scatterings 
are efficient and the attenuation due to yy annihilation is still 
not too strong. 

This parameter study aims at identifying the physical con- 
ditions in the comoving frame leading to an intense high en- 
ergy emission. Fig. [T] shows the distributions of the parameters 
corresponding to cases with an efficient inverse Compton emis- 
sion (more than 50 % of the radiated energy is due to inverse 
Compton scatterings) and a limited yy annihilation (negligible 
attenuation at the peak of the inverse Compton component). The 
most intense high-energy components are obtained for low val- 
ues of the magnetic field (B' % 100 G) and for intermediate val- 
ues of the electron minimum Lorentz factor (Fm ^ 100 - 1000), 
of the adiabatic time scale (f^x - 3 - 30 s) and of the density nf^ 
(nf^ ^ 10*^ - 10^ cm-3). 



4. Probing the parameter space of internal shocks 

In the internal shock model, the four quantities studied in 
Sect.|3] are not independent. As described in Sect.|2] they are 
determined from two sets of parameters: the first set defines 
the dynamical evolution. In the simple two shell version of the 
model, these parameters are F, k, E and t. The second set is 
related to the unknown microphysics in the shocked region: 
fB, fe, ^ and p. Therefore, we have computed 7200 spectra 
corresponding to : (i) 4 values for the mean Lorentz factor in the 
outflow, logF = 1.5, 2, 2.5 and 3; (ii) 4 values for the contrast 
which characterizes the amplitude of the variations in the initial 
distribution of the Lorentz factor in the outflow, k - 2.5, 5, 7.5 
and 10; (iii) 6 values of the injected kinetic power during the 
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relativistic ejection, log^fi/l erg s = 50, 51, 52, 53, 54 and 
55; (iv) 5 values for the variability timescale, log(T/l s) = -2, 
-1, 0, 1 and 2; (v) 3 values for the fraction of the dissipated 
energy which is injected in the magnetic field, logee = -3.5, 
-2 and -0.5; (vi) 5 values for the fraction of electrons that are 
accelerated, log^ = -4, -3, -2, -1 and 0. 

The moderate efficiency of the conversion of kinetic energy 
into internal energy by internal shocks imposes that a large 
fraction eg of this dissipated energy is injected in relativistic 
electrons to maintain a reasonable total efficiency. Therefore we 
fix fe = 1 /3. In the example presented in Fig.lT] about 7 % of the 
kinetic energy is converted in internal energy by shock waves. If 
electrons are radiating efficiently, about 2 % of the initial kinetic 
energy will be radiated. 

We also assume a slope p - 2.5 for the electron distribu- 
tion, except where mentioned otherwise. This new set of spectra 
will allow us to identify which properties of the outflow deter- 
mine the shape of the high energy spectrum, and therefore help 
to identify physical diagnostics for future Fermi data. 

4.1. The spectral shape of internal shock emission 

The effect of the six parameters , k, E, t, £B,i) on the emitted 
spectrum is now studied. We define again a "reference case" by 
f = 300, K ^ 4, E ^ 10^2 ergs-', r = 1 s, eg = 1/3 and 
^ = 10"^. Such a set of parameters corresponds to a "typical" 
GRB pulse with a peak energy Ep ^ 200 keV (source frame) due 
to the synchrotron radiation. Fig. [8] shows the evolution of the 
observed spectrum when one of the parameters is varied, while 
all other parameters are maintained constant (assuming a redshift 
z= 1). 

Maximum radius to maintain a liigli radiative efficiency. At 

very large distances from the source, the density becomes very 
low as well as the magnetic field. This increases the synchrotron 
timescale. In an equivalent way Tc is increasing and, at some 
maximum radius, can become of the order of Tm which strongly 
reduces the radiative efficiency. For such high radii, inverse 
Compton scatterings are rare due to a low density. This limit 
can then be evaluated by taking into account synchrotron radia- 
tion only. From Eqs. (|3]l, (|6]l and ( fTSl ) the condition > Fc 
leads to 

/ (/7-2)/(p-l) \ les\ ( ^ r' 

\ 1/3 ]\\I3]\\I3]\\Q-V 

Ki!)(To4^^)(3So)'(i7r^^-0^i«"' ^''^ 

with /(/<•) - {k + \)^{k - l){y/K - \)^Ik^I^. For our reference 
case, when varying one parameter only, this condition leads to 
a maximum bulk Lorentz factor F ^ 1200, a minimum contrast 
/f ^ 1 .2, a minimum injected kinetic power E ^ 1 .0 x 10"*'^ erg/s, 
a maximum timescale r S 1000 s and a minimum fraction 
of the energy injected into the magnetic field tB ~ 3 x 10"^. 
The numerical results shown in Fig. [8] agree well with these 
analytical estimates. 

Note that an additional condition should apply to 
limit the maximum radius of internal shocks (see e.g. 
iDaigne & Mochkovitchl2007h : most collisions should occur be- 



fore the deceleration radius, otherwise the propagation of the re- 
verse shock in the relativistic outflow will suppress the internal 
shock phase. For reasonable estimates of the external density, 
the deceleration radius is of the order of lO"" - 10'^ cm. From 
Eq. ([T]i, this leads to a new constraint 




with g{K) = K^l {{k -1){k+ \ 



IVIaximum density to liave an optically thin medium. On the 
other end, if internal shocks occur close to the central source, 
the density will be high. The Thomson optical depth due to the 
ambient electrons and the pairs produced by yy annihilation can 
then make the outflow optically thick. For our reference model, 
when varying one parameter only, < 1 leads to a minimum 
value of F ^ 1 10, a maximum value of /<• $ 13, a maximum value 
of £ ^ 2 X lO^"' erg s"' and a minimum value of r S 0.008 s. 

The synchrotron component at low energy. As the scaHng 
given by Eq. (l25T l for the synchrotron peak energy is quite ac- 
curate, it is not surprising to find that in mo st cases, the position 
of the synchrotron peak is simply given by (iBarraud et al.ll2005l) 

^1/2^2 ^1/2^2 ^1/2 

£syn,obs ^ T^B'Tl oc F, -^pl'^el'^ oc -2^0) (k) (43) 
with 

(V^-l)' (^-1)(^+1)3 

<i>(K) = ^ 1 . (44) 

k' 

As predicted, the observed photon energy of the synchrotron 
peak increases (see Fig.|8]l when (i) the Lorentz factor f 
decreases; (ii) the contrast k increases; (iii) the injected kinetic 
power E increases; (iv) the duration of the ejection t decreases; 
(v) the fraction increases; (vi) the fraction ^ decreases. 
This confirms that a low fraction of accelerated electrons is 
necessary to have a synchrotr on peak in the gamma-ray range 
(iDaigne & MochkovitchI 19981) and that X-ray flashes and X-ray 
rich gamma-ray bursts can be produced by internal shocks 
within "clean fireballs", i.e. outflows having a high Lorentz 
factor r (small bar yonic pollution) and a small contrast k 
(IBarraud et alj|2005h . 

There are two situations when this scaling for the syn- 
chrotron peak is not valid anymore: 

- if synchrotron radiation occurs in slow cooling regime. This 
situation would normally be rejected due to its low radiative 
efficiency. However, the synchrotron slow cooling regime can 
be compensated by efficient inverse Compton scatterings. It has 
been shown in the previous section (§ |3]l that the scaling given 
by Eq. d25T l is not accurate in this case. Even the shape of the 
synchrotron spectrum can be modified. Such cases can be found 
for instance in panel (e) of Fig.|8]for ^ 10"^; 

- if the medium is dense enough so that the synchrotron self- 
absorption frequency is above the expected synchrotron peak. 
Such highly self-absorbed cases require a large density of rel- 
ativistic electrons. As most of the spectra shown in Fig. [8] are 
computed with a low fraction f of accelerated electrons, this 
is usually not the case. In the full exploration of the parameter 
space of the internal shock model, we find that highly absorbed 
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synchrotron spectrum can be found for ^ - I. However, in this 
case the emission detected in the BATSE range corresponds to 
the inverse Compton component. This will be discussed below 

(§113). 

4.2. Spectral components in the Fermi-LAT energy range 

Conditions for intense inverse Compton emission. From the 
study made in the previous section (§[3]l, it is expected that the 
efficiency of the inverse Compton scatterings is increased by (i) 
a moderate electron minimum Lorentz factor F^, which corre- 
sponds to internal shocks with a moderate contrast k between 
the Lorentz factors of the colliding shells, and/or a large fraction 
( of accelerated electrons; (ii) a low magnetic field B', which 
corresponds to internal shocks with a high bulk Lorentz factor 
f , a moderate contrast k, a moderate injected kinetic power E, a 
large variability timescale r, and/or a low fraction of the en- 
ergy injected in the magnetic field; (iii) a low f^^, i.e. by internal 
shocks with a high bulk Lorentz factor f , a moderate contrast 
K, and/or a short variability timescale r. This is in good agree- 
ment with the results shown in Fig. [8] However, even when the 
inverse Compton emission is efficient, the corresponding spec- 
tral component is not necessarily intense, as it can be suppressed 
by photon-photon annihilation. 

Conditions for strong plioton-plioton anniliilation. As shown 
in the previous section (§ [3]), yy annihilation is important for 
large values of the optical depth t'^'^. The reason is that the yy 
annihilation and inverse Compton (in Thomson regime) cross 
sections are of the same order To investigate this effect, we have 
considered the transparency condition t!|?' < 0.1 (see § 13. It . 
It corresponds (from Eqs. ( fT9] l. ( fTOb . Q and ([B) to low bulk 
Lorentz factors f , high contrasts k, large injected kinetic power 
E, short timescales r and high fraction of accelerated electrons 
^. However, it is also favored by a high peak energy of the syn- 
chrotron component, which is also obtained for low f , high k, 
large E and short r, but low ^ and high €b . These effects are well 
observed in Fig. [8] panels (a-d) for the dynamical parameters and 
panels (e-f) for the microphysics, where it is seen in particular 
that yy annihilation is strongest for intermediate values of 

4.3. Dominant radiative process in the keV-MeV range and 
consequences at higher energy 

From our exploration of the parameter space of the inter- 
nal shock model we find, as expected from previous stud- 
ies ( Papathanassiou & Meszaros; ll996l ; iDaigne & MochkovitchI 
Il998i: iMeszaros & Reesi i2000), that there are two classes of 
spectra, depending on the radiative process responsible for the 
prompt emission in the keV-MeV range. This energy range is 
detected for instance by instruments such as BATSE, Beppo- 
SAX, HETE-2, Integral, Swift or Fermi-GSbA. These two cases 
have very different behavior in the MeV-GeV range and there- 
fore Ferm/-GBM-i-LAT observations will allow us to distinguish 
between the two possibilities. 

"Synclirotron case". The synchrotron component peaks in the 
BATSE range (keV-MeV). This case is favored in internal 
shocks as it predicts pulse shapes and spectral evolution in 
GRB lightcurves that a re in better agreement with observations 
dPaigne & Mochkovitc h 1998, 2003). The "synchrotron case" is 
found in most spectra plotted in Fig. [8] In this case, the inverse 
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Fig. 9. Effect of the slope of the accelerated electron distri- 
bution in the "inverse Compton case". The observed spectrum 
(assuming z - Y) obtained including all radiative processes is 
plotted for f = 300, /f = 4, £ = 10^^ erg s-\ r = 1 s, = 10"^ 
e^- 1/3, ^ = 1 and three different values of p. The first inverse 
Compton peak is well defined for /? > 3. 



Compton component peaks at higher energy (MeV-GeV range). 
These spectra are characterized by a high magentic field and by 
a low fraction ^ of accelerated electron^ which allows high 
values of the electron Lorentz factor r,n. In the "synchrotron 
case" most inverse Compton scatterings occur in Klein-Nishina 
regime. This leads to four types of spectra at high energy (LAT 
range): 

1 . a strong second peak with a large yy attenuation. This case 
is found for example for e^, ^ 0.03 in Fig.|8] panel (e); 

2. a weak second peak with a negligible yy attenuation. This 
case is found for example for F ~ 450 - 600 in Fig.|8] panel 
(a); 

3. a weak second peak with a strong yy annihilation. This case 
is found for example for F ~ 250 - 450 in Fig.|8] panel (a); 

4. no second peak, the high-energy emission is only the tail of 
the synchrotron component, with a cutoff in the 100 MeV - 



- This is why the "synchrotron ca se" in fast c ooling regime, which is 
our preferred case, is disfavored bv lKumar & Mc Mahon (2008). Their 
study does not consider the possibility to have ^ < 1. Therefore, the au- 
thors conclude that the "synchrotron case" in fast cooling regime is very 
unlikely, as it would involve very high contrasts k = T2/T] in internal 
shocks. With f = 1, the only possibility to reach high electron Lorentz 
factor is indeed to dissipate more energy per particle in shocks. The as- 
sumption ^ 1 made in the present study solves this problem. Note 
that Kumar & McMahon (2008) also disfavor the "synchrotron case" in 
slow cooling regime, as it implies a typical radius for internal shocks 
which is too large (of the order of the deceleration radius or larger). We 
do not discuss this case in the present study as it reduces even more 
the efficiency of the conversion of the kinetic energy of the outflow into 
radiation by internal shocks, which is already low in the fast cooling 
regime. 
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10 GeV range due to yy annihilation. This case is found for 
example for F :$ 250 in Fig.[8] panel (a). 

From case 1. to case 4., the bulk Lorentz factor F is decreas- 
ing, the contrast k is increasing, the injected kinetic power E is 
increasing, and/or the timescale r is decreasing. We discuss be- 
low the corresponding possible physical diagnostics using Fermi 
data (§l4a. 



"Inverse Compton case". The synchrotron component 
peaks at low energy and the inverse Compton component 
peaks in the BATSE range (keV-MeV). This case is usually 
called "Synchrotron Self-Compton" in the literature, and 
emerges naturally in the often considered situation where 
all electrons are accele r ated ( l^ - 1). It has bee n show n by 
iPanaitescu & MeszarosI (l2000h : IStern & PoutanenI (|2004 that 
it can reproduce the steep low-energy spectral slopes observed 
in the BATSE range. Because of the bright synchrotron com- 
ponent at low energy, possibly in the optical range, this case 
has recently been proposed to explain t he prompt emission 
of the "naked eye burst" ORB 080319b jRacusin et al.ll2008t 
iKumar & Panaites"cull2008l) and more generally o f bursts with a 
bright prompt optical emission (Panaitescu 2008). The "inverse 
Compton" case is characterized by a low magnetic field and 
a high fraction of accelerated electrons. In addition, having a 
well defined IC peak in the BATSE range requires a steep slope 
for th e electron distribution, p > 3 (see also Pe'er & Waxman 
|2004|) . While, the results obtained in the "synchrotron case" are 
very similar for all slopes p > 2, we find a large difference in the 
"inverse Compton case" between spectra computed assuming 
a slope 2 < /9 < 3 or assuming a slope p > 3. It is only for 
p > 3 that the first inverse Compton peak is well defined (the 
scatterings by electrons above Fp become negligible). This is 
illustrated in Fig.|9l 

Fig. [To] shows the same parameter study as in Fig.[8j but 
for a reference set of parameters corresponding to the "inverse 
Compton case". Note that the allowed range for each parameter 
F, K, E and r is usually more limited than in the "synchrotron 
case", especially due to the requirement of a high radiative effi- 
ciency (Mrad/Mg'^'^ > 0.5), as a lower magnetic field leads to longer 
synchrotron timescales. In the "inverse Compton case", most 
scatterings occur in Thomson regime, due to a low magnetic field 
and a low minimum electron Lorentz factor leading to Fp <K Fkn 
(see Eq. (|3TTi). The condition necessary to have the possibility of 
a second scattering in Thomson regime is Fpfp ic <K 1, i.e. 



F «F*^"'^ 



200 



B' 



100 G 



-1/5 



(45) 



For most parameters in the "inverse Compton case" this con- 
dition is fulfilled and efficient second scatterings occur, leading 
to a second inverse Compton component at high energy {Fermi 
range). The first inverse Compton component is never affected 
by yy annihilation. Therefore the spectra differ again mainly by 
their high-energy component, i.e. by the intensity of the second 
inverse Compton component. This intensity depends on whether 
most second scatterings occur in Thomson regime or are affected 
by Klein-Nishina corrections, and also on the strength of the at- 
tenuation due to yy annihilation. As long as Klein-Nishina cor- 
rections and yy absorption are not too strong at very high energy, 
the synchrotron, first and second inverse Compton components 
have relative intensities 1 : F : Y^, where Y has to be large to 



avoid that most of the energy is radiated in the synchrotron com- 
ponent in the sub-keV range. Therefore, it is difficult to avoid 
that most of the energy is radiated in the MeV-GeV range. The 
isotropic equivalent radiated energy in the BATSE range is typ- 
ically /irad.BATSE ~ 10^' - 10^"* erg. If the Compton parameter 
is F ~ 10 or more, the resulting total radiated energy is greater 
than /irad ~ 10^^ - 10^^ erg. This can lead to a crisis for the GRB 
energy budget and is another reason to di sfavor the " invers e 
Compton case" as pointed out recently by iPiran et al] ( l2008l) . 
High magnetic field can lead to smaller values of Y but most of 
the energy is radiated in the synchrotron component in this case. 
Having the first inverse Compton peak dominant requires to fine- 
tune fB ■ Moreover the peak energy of the first inverse Compton 
component has a stronger dependence on the variations of the 
physical conditions in the shocked regions (compare Eqs. ( 125] ) 
and ([30]l). Thus the "inverse Compton case" also predicts a faster 
spectral evolution during GRB pulses than in the "synchrotron 
case" and is therefore disfavored by t he observed pulse shape 
and spectral evolution in BATSE bursts dPaigne & MochkovitchI 
Il998ll2003h . 



4.4. Physical diagnostics from Fermi observations. 

As can be seen from this study, the high-energy emission 
component is shaped by several physical parameters of the 
internal shock model. It is therefore difficult to identify simple 
diagnostics that could be applied to forthcoming Fermi data. It 
is only a detailed spectral fitting covering a broad spectral range 
that will allow us to measure fundamental quantities which are 
still largely unknown for GRBs (e.g. the radius and the Lorentz 
factor of the emitting material, the typical Lorentz factor of 
radiating electrons or the magnetic field in the shocked region). 



Diagnosing tine dominant radiative process and tine pliysi- 
cal conditions in tine sliocked region. As seen in Fig.[TT] one 
can distinguish between the "synchrotron case" and the "inverse 
Compton case" from the spectral shape and then identify the 
dominant radiative process. This requires however a broad spec- 
tral range, like the one available with GBM+LAT. More precise 
informations about the physical conditions in the shocked region 
can be obtained from such observations using the following pro- 
cedure: (a) assume microphysics parameters (the initial choice 
is suggested by the general spectral shape, for instance a high 
eB and a low fraction ^ if the "synchrotron case" without bright 
IC component at high energy is favored); (b) estimate r from 
the observed lightcurve; (c) vary F and for each F adjust k and 
E to match the correct observed peak energy and fluence in the 
BATSE range; (d) determine F from the spectral shape at high 
energy. This procedure is illustrated in Figs. [TT] and [12] Good 
quality high energy observations can help in reducing the uncer- 
tainty related to the assumptions made for cb and 



IVleasuring tine Lorentz factor of tine outflow. It has been 
proposed by several authors to measure the Lorentz factor 
of the relativistic outflow from the position of the cutoff in 
the high energy spectr u m due to yy annihilation (see e.g . 
Baring & Harding 1997; Lithwick & Sari '2001; 'Baring! 12000; 
Gupta & Zhang 2008; Murase & loka 2008). It is assumed that 
three quantities have been measured: the isotropic equivalent ra- 
diated energy Ziiad.iso, the peak energy fip.obs and the typical pulse 
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duration Tobs = (1 +z)t. From Sect.lT] the radius where the emis- 
sion takes place is related to the unknown Lorentz factor of the 
emitting region r» by = 2k^/{k^ - l)rjcr ^ 2rJ ct (the differ- 
ence is less than 10% for k > 3). The comoving photon density 
can then be estimated as a function of F, : 



F.£„ 



rad,iso 



F.£' 



(46) 



where -/3 is the observed slope of the photon spectrum in the 
BATSE range and A' is the comoving width of the emitting 
region. Using the Dirac approximation for the yy annihilation 
cross section jGould & Schrede^ll967h . one gets the following 
expression for the optical depth at observed energy Eobs : 



Tyy (fobs) - CTtA 



F, [m^c^f 
' n' 

(1 H-z)£'obs 



2 \ 



E' = 



F. (mec^) 

(1 +Z)£obs 



(47) 



The cutoff energy can be estimated from the condition 
Tyy (^cuLobs) ~ 1, which givcs 



1 +z 



(1 -Hz) a-jEad ,iso 



2fii/(l +z)£p,obs 



(48) 



This estimate of the observed cutoff energy iicut.off is plotted as 
a function of the Lorentz factor F» of the emitting material in 
Fig.[T3]for the same observed quantities as in Fig.[TT] i.e. z = 1, 
firad^so = 5x10^' erg, fip.obs = 200 keV and Tobs = 2 s, and 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




assuming a slope /3 - 2.3 (note that - for most of the values 
of the cutoff energy iicut.obs plotted in Fig. [13] - the typical 

energy ~ FJ (irigC^^ /( 1 + z)^/£'cut,obs of low-energy photons that 
annihilate preferentially with photons at energy /icut.obs is above 
£p,obs, which justifies our choice of /3). On the same figure the 
measured values of iicut.obs are plotted for all spectra shown in 
Fig HU i e- for three different cases regarding the microphysics. 
This cutoff energy is measured as the energy where the slope of 
the photon spectrum falls below -3. It appears that the precise 
calculation agrees well with the approximate expression given 
by Eq. ( |48] l in the "synchrotron case" (except for a normalization 
factor of about ~ 3, which could be improved with a more 
accurate description of the low energy spectrum in equation l46b 
but that there is a large discrepancy in the "inverse Compton 
case" where Eq. (l48T l overestimates /icut.obs- The yy annihilation 
in this case is indeed negligible and the cutoff observed at high 
energy is due to the limitation of inverse Compton scatterings 
by Klein-Nishina corrections. 

From this study, it appears that using estimate of the Lorentz 
factor from the observed cutoff energy such as Eq. (|48] | is 
acceptable when the high-energy spectrum does not show any 
, new bright component in addition to the low-energy spectrum. 
To confirm any value of the Lorentz factor measured by this 
method, one should do a detailed modelling of the spectrum 
to check whether the process responsible for the observed 
cutoff has been correctly identified. Note that the emission 
from electron-positron pairs created by yy annihilation is not 
included here and tha t it could provide additional diagnostics 
urase & lokal l2008l ). Note also that the final shape of the 
spectrum above the cutoff energy can become complicated when 
considering the spectrum integrated over the whole duration of 
a pulse, as we show in the next section (see also the detailed 
discussion by iGranot et al.l 120081) . It can make the diagnostics 
more difficult and emphasizes again the necessity of a detailed 
modelling. 



5. Observed time profiles and spectra 

All the spectra shown in the previous section are computed for 
a "typical collision", using the two shell version of the internal 
shock model. However even for a single pulse, the emission ra- 
dius, and the physical conditions in the shocked medium, can 
span several orders of magnitude during the propagation of the 
"internal shock" waves (see Fig.[T]). This leads to a spectral evo- 
lution during the observed pulse that is entirely missed by the 
two shell model. In the present section, we show examples of 
synthetic bursts computed by coupling the detailed model for the 
dynamics of the relativistic outflow (§ Sect. 12. 11 ) with the radia- 
tive code (§ Sect. |2.3l ). and we discuss the predicted spectral evo- 
lution, as well as the high energy emission (Ferm/-LAT range). 
In all examples, a redshift z = 1 is assumed. 



Fig. 13. Measuring the Lorentz factor of the emitting region 
from the observed cutoff energy. The cutoff energy /Scut.obs due 

to yy annihilation is plotted as a function of the Lorentz factor 
F, of the emitting region, as predicted from Eq. ( l48T l, assuming 
z = 1, -Erad^so = 5 X 10^' erg, Ep^obs = 200 keV and /3 = 2.3. The 
cutoff energies measured in all the spectra shown in Fig. [TT] are 
indicated with four different symbols corresponding to the four 
different assumptions regarding the microphysics parameters. 



5.1. A single pulse burst 

We present three synthetic single pulse bursts corresponding to 
the same relativistic outflow : a total duration of the relativistic 
ejection phase = 2 s, a Lorentz factor varying from 100 to 
400 during the ejection as in Fig. [1] and E - 5 x 10^^ erg s '. 
The dynamics is computed using a discretization of the outflow 
in 1000 shells. The three cases differ by different sets of micro- 
physics parameters : (i) "synchrotron case" with a high magnetic 
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field (see Fig.IHi, eg = ee = 1/3, f = 3 x 10"^ and = 2.5 ; 
(ii) "synchrotron case" with a low magnetic field (see Fig.fTSll. 
fB = 5 X 10"\ ee = 1/3, f = 2 X 10"^ and p = 2.5 ; (iii) "inverse 
Compton case" (see Fig.[T6]l. eb = 10"^, ee = 1/3, ( - 1 and 
p - 3.5. In this last case, the contrast k - Fmax/rmm of the 
initial distribution of the Lorentz factor has been increased (F 
varies from 100 to 600 instead of 400) to increase the dynamical 
efficiency, and thus compensate for a lower radiative efficiency 
(as well as a lower fraction of the emission that is radiated 
in the Fenni-GBM range). The three pulses have comparable 
isotropic equivalent energies radiated in the GBM range. The 
lightcurves in the GBM+LAT range are plotted for each case, 
as well as the time-integrated spectrum during the rise, the 
decay and the whole duration of the pulse. GRB lightcurves in 
the keV-MeV range usually show a hard-to-soft evolutio n (see 
e.g. iBhat et al1fT994t lFbrdetaDfT99l iNorris et al.lll996h . This 
spectral evolution is found in these three examples of synthetic 
GRBs, as the spectrum during the rise peaks at higher energy 
than during the decay phase. 

In the GBM range, the three lightcurves are quite similar, 
except for a faster pulse decay in the inverse Compton case. 
We checked in th e three cases tha t the expected spectral 
evolution (see e.g. iNorris et"aL I [19961) in the keV-MeV range 
is r eproduce d, in agreement with the previous results of 
iDaign e & MochkovitchI (Il998l) : the photon flux peaks earlier at 
higher energy and the duration of the pulse increases at lower 
energies. 

In the LAT range on the other hand, the spectral evolution 
and the corresponding behavior at high energy are different in the 
three considered cases. In the "synchrotron case", the physical 
process responsible for the radiation is not the same in the GBM 
range (synchrotron) and in the LAT range (synchrotron+possible 
additional inverse Compton component, depending on the in- 
tensity of the magnetic field). Therefore, the lightcurves at low 
and high-energy do not look similar. In particular, the inverse 
Compton component at high energy emerges later than the syn- 
chrotron component, increasing the duration of the pulse in the 
LAT range. It is due to an evolving Compton parameter during 
the pulse duration. This effect is more important when inverse 
Compton scatterings become more efficient (compare Figs. [T4l 
and [15]). For the lowest values of eB, the lightcurve at high en- 
ergy could even peak with a delay with respect to the lightcurve 
in the GBM range if the inverse Compton component becomes 
more intense that the synchrotron component in the high energy 
range. This behavior of the LAT lightcurves is due to the evolu- 
tion of the physical condition in the shocked medium along the 
propagation of the shock wave. It is illustrated in Fig. [IT] where 
the ratio of the inverse Compton over the synchrotron compo- 
nent is plotted as a function of the observer time, as well as the 
dynamical timescale t'^^ and the synchrotron timescale t'^^^ : 

- Initially, due to the shape adopted for the initial distribution 
of the Lorentz factor in the outflow, the shock is weak and the 
dissipated energy per particle is low. This results in moderate 
electron Lorentz factors Fm, and therefore large synchrotron 
timescales. On the other hand, these early times correspond 
to small radii so the dynamical timescale is still small. In this 
first phase, f^^^ i f^x ^i^d efficiency of inverse Compton 
scatterings is large, as a large fraction of the shocked region 
is populated by relativistic electrons (see Sect.[3j. It results 
in a weak precursor in the GeV lightcurve. This precursor 
can disappear if a different initial distribution of the Lorentz 



factor in the outflow is adopted, especially if it leads to an 
immediate violent shock (for instance with an initial discon- 
tinuity). 

In a second phase (around the peak of the pulse in the GBM 
range), the shock becomes stronger, Fm increases and the 
synchrotron timescale decreases. On the other hand, as the 
radius increases, the dynamical timescale increases. This re- 
sults in fsy„ <s: and a low efficiency for inverse Compton 
scatterings. The emission at high energy is dominated by the 
synchrotron component. 

A late times (tail of the pulse), the synchrotron timescale 
increases again (mainly due to the decrease of the mag- 
netic field as the radius increases) and the efficiency of 
inverse Compton scatterings increases again. The inverse 
Compton component becomes dominant again in the GeV 
range, which results in a prominent GeV tail of the pulse. 



On the other hand, in the "inverse Compton case", inverse 
Compton scatterings are the dominant radiative process both 
in the GBM and the LAT range. Therefore, the lightcurves are 
much more similar in the different energy bands. 

It appears clearly from Figs. [141 [15] and [16] that in addition 
to the time-integrated spectrum, the observed spectral evolution 
and the comparison of the GBM and LAT lightcurves are also 
powerful tools to diagnose the dominant radiative process and 
the physical conditions in the shocked medium (electron distri- 
bution and magnetic field). 




Fig. 17. A single pulse burst in the "synchrotron case" with 
a low magnetic field: dynamical evolution during the pulse. 

The evolution of the synchrotron timescale, the dynamical 
timescale and the ratio of the energy radiated in the inverse 
Compton over the energy radiated in the synchrotron component 
is plotted as a function of the observer time for the synthetic sin- 
gle pulse burst shown in Fig. [15] 
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Fig. 18. An example of a multi-pulses burst : initial distribu- 
tion of the Lorentz factor in the relativistic outflow. The initial 
Lorentz factor in the outflow is plotted as a function of the ejec- 
tion time (top axis) or equivalently the distance from the source 
at the end of the ejection phase (bottom axis). 



5.2. A multi -pulses GRB 

Fig. [19] shows two examples of more complex synthetic GRBs. 
The dynamical evolution is the same in both cases, assuming the 
initial distribution of the Lorentz factor plotted in Fig. [18] which 
leads to 4 main pulses in the lightcurve. It is computed using 
a discretization of the outflow in 4000 shells. The two exam- 
ples correspond to two different sets of microphysics parameters 
("synchrotron case" with a high or a low magnetic field). These 
examples illustrate that in a complex burst, each pulse exhibits a 
hard to soft evolution in the main spectral component and that - 
when possible - the spectral analysis should be made by integrat- 
ing the spectrum over a pulse rather than over the whole duration 
of the GRB. The spectral evolution and the behaviour at high en- 
ergy that were identified for single pulse burst are also observed 
in these multi-pulse GRBs. In particular, the lightcurve above 1 
GeV in the "synchrotron case" with a low magnetic field shows 
more flat topped pulses and prolonged emission in the pulse de- 
cays. 

6. Conclusions 

We have developed a detailed model for the prompt emission of 
gamma-ray bursts in the framework of the internal shock model. 
It combines dynamical simulations that follow the evolution of 
the physical conditions (Lorentz factor, density and energy den- 
sity) in the shocked regions along the propagation of "internal" 
shock waves in the relativistic outflow, and a time-dependant 
radiative code to compute the emission from shock-accelerated 
electrons, including the most relevant processes (adiabatic 
cooling, synchrotron radiation and self-absorption, inverse 
Compton scatterings and photon-photon annihilation). We have 
used our model to explore the parameter space of the internal 
shock model and identify physical diagnostics for Fermi data. 



We list here our main results: 

1 . The comparison of the results of the radiative code with stan- 
dard analytical or semi-analytical estimates of the emitted 
spectrum shows that the synchrotron component is usually 
well predicted by the analytical spectrum from ISari et all 
(Il998l) . except when inverse Compton scatterings become 
efficient. 

2. We show that the prediction of the high-energy com- 
ponent based on the time-averaged electron distribution 
and the time^averaged synchrotron photon spectrum (e.g. 
ISari&EsiDll2()0lh is less precise, especially above the peak 
of the inverse Compton component. The accuracy of this 
time-averaged prediction decreases as the efficiency of in- 
verse Compton scatterings increases. This is mainly due to 
time-dependent effects : time-averaged distributions do not 
take into account the time needed to build the photon field. 

3. An important consequence of the previous result is that 
the electron cooling rate, and therefore the time-averaged 
electron distribution, is not cor rectly pred icted by the stan- 
dard analytical estimate from IS ari et al] fl998) when in- 
verse Compton scatterings become important. This affects 
the spectral shape of the synchrotron component as well. 
We find that the resulting spectral slope (photon spec- 
trum) below the peak energy is steeper than the usual value 
-3/2 in synchrotron f ast cooling regime, in agreement with 
iDerishev et alj ( l200lh who shows that the slope can be- 
come as steep as -1 when Klein-Nishina corrections are 
important. This may reconcile the synchrotron radiation 
with the o bserved d istributio n of the low-en ergy slope a in 
BATSE (.Preece et al...2000; Kaneko et al.l2006) and HETE-2 
dSakamoto et alj|2005h bursts. We will investigate this ques- 
tion in a forthcoming paper 

4. When exploring the parameter space of internal shocks, we 
find that two classes of broad-band spectra can be expected, 
which correspond to different physical conditions in the 
shocked region : 

- "Synchrotron case", where the dominant process in the 
Fermi-GBM range is synchrotron radiation. It requires 
high electron Lorentz factors and therefore implies that 
only a fraction of the electrons are shock-accelerated. 
The intensity of the inverse Compton component in the 
LAT range depends on the intensity of the magnetic field 
but remains always limited due to Klein-Nishina cor- 
rections. A high-energy cutoff is present due to photon- 
photon annihilation, which allows us to estimate the 
Lorentz factor of the emitting material. As the high- 
energy lightcurve is made of two different components 
(synchrotron -i- inverse Compton) whose ratio evolves 
during a pulse, the lightcurves can differ in the GBM and 
the LAT range. In particular, for low fB/fe ratios that fa- 
vor inverse Compton scatterings, a pulse in the GBM can 
be followed by a tail of GeV photon in the LAT. When 
there is a significant inverse Compton component in the 
LAT, the pulse above 1 GeV can even peak with a delay 
compared to the GBM. 

- "Inverse Compton case", where the synchrotron compo- 
nent peaks at low energy and the dominant process in 
the GBM range is inverse Compton. A second inverse 
Compton peak is usually present at higher energy but 
further scatterings are suppressed by Klein-Nishina cor- 
rections. The cutoff at high energy can be either due to 
yy annihilation or to the Klein-Nishina suppression of 
inverse Compton scatterings, which makes more difficult 
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Table 1. Physical diagnostics for Fermi data: 



Observational constraints : 








-Is there an additional component in the LAT spec- 


no (or very weak) 


yes 


yes (possibly very intense) 


trum ? 








-Is the LAT lightcurve in GeV range prolonged or 


yes (tail) 


yes (delayed peak) 


no (or very short delay) 


delayed compared to GBM lightcurve ? 








Consequences: 








Radiative processes: 








Dominant process in GBM range 


synchrotron 


synchrotron 


IC (first peak) 


Dominant process in LAT range 


synchrotron + very weak IC 


synchrotron + IC 


IC (second peak) 


Conditions in shocked regions: 








Electron distribution 


high Tm, low ^ 


high Tni, low ^ 


moderate Tn,, — > 1 


Magnetic field 


high (6b ^ 1/3) 


low (6b ^ 1) 


low (6b <c 1) 


Properties of the relativistic outflow: 








Can the Lorentz factor in the outflow be measured 


yes 


yes 


? (see text) 


from the high-energy cutoff (assuming it is due to 77 








annihilation) ? 









These diagnostics are based on the results shown in Sect.|4]and|5] 



to estimate the Lorentz factor of the emitting material 
in this case. A steep slope for the electron distribution 
is necessary to have two well distinct peaks in the in- 
verse Compton component. The relative intensity of the 
synchrotron and the two inverse Compton peaks in the 
spectrum depends on the ratio e^je^. It requires some 
fine-tuning to have a dominant first inverse Compton 

peak in the GBM range, which can lead to a energy cri- 

f 1 1 ' 

sis (Piranetal. 2008). As the emission detected in the 

GBM and LAT lightcurves is due to the same process (in- 
verse Compton), the lightcurves in the two instruments 
are much more similar 
5. This study allows us to define physical diagnostics for Fermi 
data, based on the spectral shape and the spectral evolution, 
that are summarized in Tab.[T] We plan to apply these diag- 
nostics as soon as Fermi GRB data will be made public. Our 
study emphasizes however that a detailed broad-band spec- 
tral modelling is always necessary to reach firm conclusions 
regarding the properties of the outflow and the physical con- 
ditions in the shocked regions. 

Several arguments (energetics, spectral evolution) already favor 
the "synchrotron case". Bursts detected by Fermi both by the 
GBM and the LAT should allow us to firmly distinguish these 
two possibilities and show if all bursts are in the "synchrotron 
case" or if the "inverse Compton case" can also be found. 

In the "synchrotron case", it is assumed that only a frac- 
tion of the electrons is accelerated to very high Lorentz factors 
(Fni ~ 10"^). If a small fraction of the dissipated energy in the 
shock is injected in the remaining electrons, they will have a 
Maxwellian distribution with a mean Lorentz factor of a few. We 
will investigate in a future work what could be the contribution 
of these electrons to the emission and more generally discuss 
the prompt optical emission of GRBs in the framework of the 
internal shock model. 



Appendix A: Radiative processes 

We list here the source and loss terms appearing in the equations governing 
the evolution of the electron distribution (Eq. )15t ) and of the photon spectrum 
(Eq.Ol). 



A.1. Adiabatic cooling 

dy I y 
dt' lad tpx 



A.2. Synchrotron emission 



(A.1) 



We assume an isotropic distribution for the pitch angle a between the electron 
velocity and the magnetic field, so that ^sin~ = 2/3. The synchrotron power 
is given bv lRvbicki & LightmanI i 19791) and leads to 



dy I 

Isyn 



B y 



3e 



B<I> 



Ay Ay)) ' 



where the synchrotron frequency is defined by 

Vsyn (r) = 



3e 7 



and the function <D is defined by 

9V3 r+" 

This definition has been chosen so that dx <D(j:) = 1 . 

A.3. Synclirotron seif -absorption 

The cross-section is given bv lRvbicki & LightmanI 09791) : 

1 nr(r)| 



(A.2) 
(A.3) 

(A.4) 

(A.5) 



o-sa (r, v') -- 



d\nn , ,^ 
2-——(y,t') 
amy 



(A.6) 



In the present version of our radiative code, the corresponding source term 
dyldt'\^^ is not included in the equation for the evolution of electrons (Eq. (Ts))- 
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A.4. Inverse Compton scatterings 

We use the kernel derived bv lJonesHl968l) , which has an excellent accuracy, even 
in the Klein-Nishina regime. 

= -T — - dv'y' —n-At')K(yy,v') , (A.7) 
dt' lie 4 meC 72 J J v' 
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with 



— hcTTC 

4 y 



v' fdv' 

yl J v' 



)K{'y,v',v') 



(A.8) 



K('y,v',v') = !-r if < e < ? , 



2q\nq + (l+2q)(i-q)+^{l-q)p^ 
2 1 + 4yeq 



if e < e < 



1 + 4y£ 



(A.9) 



where e = hv'/m^c^, e = hv' jm^c- and q = e/4ye/(y- e) . In the present 
version of the code, the corresponding loss term 



Aly' I 

dt' lic.loss 



'('') ^ dy«(y, t')a;^(y,v') 



is not included in the right-hand part of the equation for the evolution of photons 
(Eq. (16)). For this reason, we can not compute the emitted spectrum when T^'^ S 
1 and can not reproduce a Comptonized spectrum. 



A.5. Photon-photon annihilation 



We use the exact cross section given by Gould & Schredei'^1967D for an isotropic 
photon field, which is a good approximation as long as the radiative timescale is 
small compared to the dynamical timescale (see Sect. |2.3K 



CTyy(v',v') 



where 



1+/J2 



l-/?2 



-2/5 



1+/; 



<I+/3)/(l-« 

— ln(l , 



(A. 10) 



hv'hv' - {msC^y^ 
hv'hv' 



In the present version of the code, the source and loss terms due to pair pro- 
duction and pair annihilation are not included in Eq. il5\ . and the associated 
radiation is not included in Eq. (16). 



Appendix B: Numerical method 

To solve the system of the two Eqs. (15) and )16t for the evolution of the electron 
distribution n (y, ;') and the photon spectrum nyi(t') in the comoving frame of the 
shocked material, we have developed a numerical scheme that is a good com- 
promise between accuracy and computing speed. We use normalized variables 
t = t' /t'^^, e = hv'/nicC-, h(y,t) = n(y,i) jnf^ 3nAh^{t) = (nicC^ I li^ iiyi (t' ) I rif^ . 
The photon spectrum at time after the k"" step is stored in a fixed grid ej for 
7 = 1 — ► M, i.e. = «irj (fit). At F= 0, no photons are present so that nj(0) = 0. 
To take advantage of the short radiative timescale, the electron distribution at 



time fi; is stored in a moving ("Lagrangian") grid y^^^ = y\ (fjj) for i = 



An. 



dy n{y,i) . 



This insures that the number of electrons is exactly conserved. At F = 0. this 
grid is adjusted to the distribution of accelerated electrons given by Eq. (5)' 



(0) 



I m. r 



,(0 



Fm and Ah 



,(0 



\t-p 



In practice, all 

,(0) 



results presented in this paper are obtained using logaiithmic grids for ej and y: 
with N = M = 100. The time-evolution of nj''* is computed using a discretized 

version of Eq. (16) . On the other hand, to compute the time-evolution of y*', we 
replace Eq. (Ts) by 



' I 

, df lad 



dy I 



dy I 



(B.l) 



In fast cooling regime the values of y|''* decrease rapidly, which is a valuable ad- 
vantage, as the timestep of the simulation is given by the radiative timescale of 
the most energetic electrons still present in the medium. We observe a rapid in- 
crease of the timestep following the cooling of the electrons, and reach therefore 
f = 1 in a small number of steps. In addition, as -dy/dtis an increasing function 
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Fig. C.l. Time evolution of the Compton parameter. The 

Compton parameter Y(t') is plotted for r,n - 100, B' - IC* G, 
fg^ = 10 s and four values of the initial density of relativis- 
tic electrons which are labelled in the Figure. The numerical 
calculation (thick solid line) is performed including adiabatic 
cooling, synchrotron radiation (without self-absorption) and in- 
verse Compton scatterings. The analytical expression of F(f') de- 
scribed in appendix|C]is plotted with a dashed line. A horizontal 
thin line indicates the limit Y - \ between a regime where elec- 
tron cooling is dominated by synchrotron radiation or inverse 
Compton scatterings. A vertical dotted line indicates the syn- 
chrotron timescale for electrons with Lorentz factor Fm. 



of y, the width of each cell in the electron grid (y*' - y*'| ) is decreasing with 
time, which allows at the end of the simulation to merge cells of cooled electrons 
and therefore gain an additional reduction of the computing time. Energy is not 
necessarily conserved in our code due to the finite size of the electron and photon 
grids: typically, the en'or is a few percents and is never lai'ger than 10 %. 



Appendix C: Formation of the radiation field: the 
Compton parameter Y(t'). 

When only adiabatic cooling and synchi'otron radiation are considered, the equa- 
tion for the evolution of electrons has a simple analytic solution, given by 



n(y,0 = [^:^^] „(TO(y,O,''=0) 



(C.l) 



where yo (y, f') is the initial Lorentz factor at f' = of an electron having cooled 
to the Lorentz factor y at time f' , i.e. 



ro (y, t') ■- 



Ec 



At time f, electrons have Lorentz factors in the interval 

Ec „ ,„ Fe 



r,„(f') 



1 + . 



< y < ruW) ' 



^^^e''/'» + l 
The time-averaged electron distribution 

7" r"*'H(y,r') 

'cx Jo 



(C.2) 



(C.3) 
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can be computed exactly from these expressions. An accurate approximate ex- 
pression is given by Eq. i22\ . When inverse Compton scatterings are included, 
this exact solution for the electron evolution is not valid any more, but remains 
very accurate as long as the synchrotron process is still dominant, i.e. as long as 
Y(t') <K 1. It is therefore interesting to compute the evolution of the Compton 
parameter in this case. It is defined as F(f') = Hrad('')/ ( B' " /Sjr). After some 
algebra, it can be written as 



y(t') = T(p-i)('^T«rc?sy„(r„))r^ 



^ m 



'Jo 



dyy" 



y/4 



(1 +).)(//'« -i) 



(C.4) 



The final value at Fex = Y{t'= 4) is given in Sect.[32] We plot in Fig. lcTl the 
evolution of Y(t') for increasing values of Fe^. The agreement between the nu- 
merical result and the analytical solution is excellent as long as the synchrotron 
radiation is still the dominant process, i.e. for Y(t') < 1. The discrepancy is 
larger when the transition Y' = I occurs in a time shorter than the synchrotron 
timescale. In any case, in fast cooling regime the energy density contained in the 
photon field at f = equals njl'^'^ (initial energy density injected in relativistic 
electrons) and therefore the asymptotic value is always equal to the theoretical 
prediction Y^„ = uf^ I {^B'^ if the Thomson regime is valid (which is the 
case in all examples plotted in Fig. lC.U . 
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Fig. 8. The effect of internal shock parameters on the emitted spectrum ("synchrotron case"). We use the simple two shell 
version of the internal shock model (see text) and define a "reference case" by F = 300, k = A, E = erg s"\ r = Is, 

fe = = 1 /3, <r = 10"^ and p - 2.5. Each panel shows the evolution of the observed spectrum (assuming z = 1) when one parameter 
is varied, while all other parameters are maintained constant. Two effects can limit the parameter range: electrons become radiative ly 
inefficient ("efficiency limit") or the medium becomes optically thick due to the intense production of pairs ("transparency limit"). 
For each limiting case corresponding to the transparericy Umit (panels a,b,c,d), the spectrum that would be observed without yy 
annihilation is plotted with a dashed line, (a) effect of F. The transparency limit is reached for F < 109 and the efficiency limit for 
F > 1400 ; (b) effect of k. The efficiency limit is reached for k < 1.5 and the transparency limit for k > 12.5; (c) effect of E. The 
transparency limit is reached for E > 2x 10^* erg s"' ; (d) effect of t. The transparency limit is reached for t < 0.008 s; (e) effect of 
Cfl. The transparency limit is never reached. The spectrum that would be observed without yy annihilation is plotted with a dashed 
line for = 10"^ and = 1/3; (f) effect of The transparency limit is never reached. The spectrum that would be observed 
without yy annihilation is plotted with a dashed line for ^ = 10"^. 
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Fig. 10. The effect of internal shock parameters on the emitted spectrum ("inverse Compton case")- Same as in Fig.|8]with a 
new "reference case" defined by f = 300, k = 4, E - 10^^ erg s"', r = 1 s, eB = lO^^*, fe = 1/3, ^ = 1 and p = 3.5. (a) effect of f . 
The transparency limit is reached for f ^ 120 and the efficiency limit for f ^ 290 ; (b) effect of k. The efficiency limit is reached 
for K ^ 4 and the transparency limit for k ^ 12; (c) effect of E. The efficiency limit is reached for £ ^ 1.1 x 10^^ erg s"' and the 
transparency limit is reached for E ^ 1.1 x 10^^ erg s"'; (d) effect of t. The transparency limit is reached for t ^ 0.01 s and the 
efficiency limit is reached for r ^ 0.9 s. In each panel, the observed spectrum obtained without including -yy annihilation is plotted 
with a dashed line for the value of the parameter corresponding to the transparency limit. 
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Fig. 11. Physical diagnostics from Fermi observations: spectral shape. A GRB pulse at z = 1 characterized by a peak energy 
^p.obs = 200 keV, an isotropic radiated energy in the BATSE range £'25keV-iMeV = 5 x 10^^ erg and a duration (1 + z)t = 2 s 
is modelled using the two-shells version of the internal shock model and different assumptions for the microphysics. Top-left: 
"synchrotron case" with a high magnetic field (eg - = 1/3, ( - 10"^ and p = 2.5); Top-right: "synchrotron case" with a low 
magnetic field (e^ = lO"^, eg = 1/3, ^ = 10"^ and p - 2.5); Bottom-left: " inverse Compton case" with a high magnetic field 
(efi = = 1/3, ^ = 1 and p = 3.5); Bottom-right: " inverse Compton case" with a low magnetic field (cb = 10"^, = 1/3, ^ = 1 
and p = 3.5). In each case, the evolution of the observed spectrum is plotted for an increasing bulk Lorentz factor T (calculation 
including all radiative processes), while the other parameters (k, E, t) are adjusted to match the imposed observed quantities. Spectra 
in dotted lines correspond to cases where pair creation should not be neglected (t!J?' > 0.1). 
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Fig. 12. Physical diagnostics from Fermi observations: properties of the relativistic outflow. For each case considered in Fig.fTTI 
we plotted the values of the injected kinetic power E and the contrast k that have been adjusted for each bulk Lorentz factor F to 
match the imposed observed quantities (redshift z = 1, peak energy fip.obs = 200 keV, isotropic radiated energy £25 kev-i Mev = 
5 X 10^' erg and duration (1 + z)t = 2 s). This adjustement is done within 10% uncertainty, which explains why the curves are not 
perfectly smooth. The corresponding typical internal shock radius Ri^ is also shown. The dotted lines correspond to cases where pair 
creation should not be neglected (t!^' > 0.1). 




Fig. 14. A single pulse burst in the "synchrotron case" with a high magnetic field. The dynamics is the same as in Fig.[T]except 
for E = 5x10^-' erg s"' . The assumed microphysics parameters are ee = £e - 1 /3, f = 3 x lO"-' and p = 2.5. All radiative processes 
are included in the calculation. Left: observed lightcurves in Fermi-GBM and LAT range. The synchrotron (thin solid line) and 
inverse Compton (thin dashed line) components are also shown. Right: observed time-integrated spectrum during the rise, the early 
decay and the whole duration of the pulse. 
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Fig. 15. A single pulse burst in the "synchrotron case" with a low magnetic field. Same as in Fig.[T4]except for the microphysics 
parameters: - 5 x 10'^, eg = 1/3, ^ = 2 x 10"-' and p - 2.5. 




Fig. 16. A single pulse burst in the "inverse Compton case". Same as in Fig.[T4]except for the initial distribution of the Lorentz 
factor that varies from 100 to 600 and for the microphysics parameters: €b = 10"^, - 1/3, ^ = 1 and p = 3.5. In the left panel, 
the lightcurves observed both in the Fermi-GBM+LAT energy range are entirely dominated by inverse Compton emission. 
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Fig. 19. An example of a multi-pulses burst : lightcurves. The dynamics is computed for the initial distribution of the Lorentz 
factor plotted in Fig. [18] assuming a total duration of the relativistic ejection phase = 8 s and an injected kinetic power E - 
5 X 10^^ erg s"' . The two examples differ by the assumptions for the microphysics parameters. Le/f; "synchrotron case" with a high 
magnetic field. The microphysics parameters are eB = ee = 1/3, ^ = 0.002 and p - 2.5. The observed peak energy of the spectrum 
integrated over the whole duration of the burst is ~ 180 keV. Right: "synchrotron case" with a low magnetic field. The microphysics 
parameters are ee = 8 x lO""*, eg = 1 /3, ^ = 0.001 and p - 2.5. The observed peak energy is ~ 250 keV. In both cases the lightcurves 
are plotted in different energy channels (ferwi-GBM+LAT) assuming a redshift z = 1 . 



